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THEORY  OF  OPTIMUM  SHAPES  IN  FREE -SURFACE  FLOWS 


PART  I 

Variational  Calculus  Involving  Singular  Integral  Equations 

by 

T.  Yao-tsu  Wu  and  Arthur  K.  Whitney 

California  Institute  of  Technology 
Pasadena.  California 

The  general  problem  of  optimum  hydromechanical-shapes  arising  in  a  wide 
class  of  free  surface  flows  can  be  characterized  mathematically  as  equi¬ 
valent  to  the  extremization  of  a  functional  of  the  form 

1 

n 

l[u]  =  j  £(u{x),  v(x),  x)dx 
-1 

where  f  is  an  arbitrary  function  which  is,  in  general,  nonlinear  in  u,  v, 
x,  and  the  unknown  argument  functions  u(x)  and  v(x)  are  related  by  the 
singular  Cauchy  integral 

'W  =  fr  <H<1>  • 

-1 

Application  of  the  variational  method  to  l[  u]  yie’ds,  in  analogy 
to  the  Euler  differential  equation  in  the  classical  theory,  the  following 
singular  integral  equation  of  the  Cauchy  type 

,  £  f  (u(t),v(t),t) 

fu(u(x),  v(x),x)  =  -  dt  ( |  X I  <  1 )  , 

-1 

where  f^  s  3f/3u  and  f  =  3f/3v.  This  equation,  which  is  a  necessary 
condition  for  the  existence  of  an  extremal  l[u],  combines  with  the  integral 
definition  of  v(x)  to  give  a  pair  of  singular  integral  equations,  which  are  to 


2 


be  solved  for  a,  v  under  appropriate  end  conditions  and,  possibly  under 
additional  isoperimetric  constraints.  Consideration  of  the  second  varia¬ 
tion  of  I  leads  to  the  inequality 

fuu*U’V,X^+fv^U,V,X^°  ( I x I  <  1 ) 
as  a  necessary  condition  for  the  extremum  of  I  to  be  a  minimum. 

Analytical  solutions  by  the  method  of  singular  integral  equations 
and  some  approximate  methods  by  Fourier  series  expansions  are  discussed 
for  the  linearized  theory.  The  general  features  of  these  solutions  are 
demonstrated  by  numerical  examples. 


1 .  Introduction 

In  recent  studies  undertaken  by  the  authors  it  has  been  observed 
that  the  determination  of  the  optimum  hydromechanical  shape  of  a  body  in 
a  free  surface  flow  invariably  results  in  a  new  class  of  variational  prob¬ 
lems,  in  which  the  unknown  functions  are  related,  not  by  differential 
equations  as  in  the  classical  calculus  of  variations, but  by  a  singular 
integral  equation  of  the  Cauchy  type.  These  recent  studies  include  the 
following  problems:  (i  )  the  optimum  chape  of  a  plate  planing  on  a  water 
surface;  (ii)  the  body  profile  of  minimum  pressure  drag  in  symmetric 
cavity  flows;  and,(iii)  the  cavitating  hydrofoil  having  a  maximum  lift-drag 
ratio  under  a  set  of  isoperimetric  constraints.  Indeed,  the  physical  prob 
lems  belonging  to  this  class  embrace  a  wide  range  of  interest.  Closely 
related  examples  are  numerous,  such  as:  (iv)  the  optimum  circulation 
distribution  of  a  lifting  line  in  the  aerodynamic  wing  theory;  (v)  the  hull 
shape  of  a  thin  ship  which  has  minimum  wave  resistance;  and  ,  (vi)  a  class 
of  mixed-type  boundary  problems  pertaining,  butnot  limited  to  elliptic 
paxtial  differential  equations.  Only  a  few  special  cases  from  this  general 
class  of  problems  have  been  solved,  the  optimum  lifting  line  being  an  out¬ 
standing  example. 

There  are  several  essential  differences  between  the  classical 
theory  and  this  new  class  of  variational  problems.  First  of  all,  the  "Euler 
equation"  which  results  from  the  consideration  of  the  first  variation  of  the 
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functional  in  this  new  class  is  a  singular  integral  equation  of  the  Cauchy 
type  which  is,  in  general,  nonlinear.  This  is  in  sharp  contrast  to  the 
Euler  differential  equation  in  the  classical  theory.  Another  noteworthy 
feature  of  this  new  class  of  variational  problems  is  that  while  regular  be¬ 
havior  of  solution  at  the  limits  of  the  integral  equation  may  be  necessary 
on  physical  grounds, the  mathematical  conditions  which  insure  such  behavior 
generally  involve  functional  equations  which  are  difficult,  and  sometimes 
simply  impossible,  to  satisfy. 

Because  of  these  difficulties  and  the  fact  that  only  very  limited 
techniques  are  known  for  solving  nonlinear  singular  integral  equations,  a 
theory  for  this  new  type  of  variational  problems  has  not  been  fully  developed. 
Attempts  are  made  here  to  present  some  preliminary  results  of  this  study. 
After  the  problem  is  stated  in  its  general  form,  the  variational  methods 
are  applied  in  Section  3  to  derive  two  necessary  conditions  of  optimality  of 
the  functional  from  the  considerations  of  the  first  and  second  variations. 

They  correspond,  respectively,  to  the  Euler  differential  equation  and  the 
Legendre  condition  in  the  classical  theory. 

Following  the  general  formulation,  solutions  of  the  singular  integral 
equations  are  sought  for  the  linear  case  when  the  functional  is  a  quadratic 
form  in  its  argument  functions  -  -  a  case  which  seems  to  be  the  least  dif¬ 
ficult,  and  is  expected  to  retain  the  important  features  of  the  corresponding 
nonlinear  problems.  It  is  of  interest  to  note  that  the  linear  problem  can  al¬ 
ways  be  converted  to  a  Fredholm  integral  equation  of  the  second  kind,  for 
which  a  well -developed  theory  is  available.  Moreover,  when  the  coefficients 
of  the  linear  integral  equations  satisfy  certain  relationships,  analytical 
solutions  have  been  obtained  in  closed  form,  by  the  method  of  singular 
integral  equations.  To  aid  practical  applications,  an  approximate  method 
(which  is  essentially  the  Rayleigh -Ritz  method)  employing  a  discretized 
Fourier  series  representation  of  the  desired  solution  is  discussed,  and  the 
results  are  compared  with  some  known  exact  solutions. 

2.  Statement  of  the  Problem 

The  general  optimum  problem  considered  here  may  be  stated  as 
follows:  To  find  the  real, extremal  function  u(x)  of  a  real  variable  x. 


required  to  be  Holder  continuous^  in  the  region  -1  <  x  <  I,  together  with 
its  finite  Hilbert  transform 

^  1 

v(x)=i£  (-1<x<1)  ,  (1) 

-1 

where  the  integral  with  symbol  C  signifies  its  Cauchy  principal  value,  so 
that  u  and  v  minimize  the  functional 

1 

J[u]  fo(u(x),  v(x),  x)dx  ,  (2) 

-1 

satisfy  M  isoperimetric  constraints 
1 

J|[u]  =y  fj(u(x),  v(x),  x)dx  =  const.  =  (1=1,2,  .  .  .  ,  M) 

_1  (3) 

and  satisfy  the  conditions  near  the  end  points  x  =  ±  1, 

or, +ip, 

u(x)  =  u^(x)/(x-ck)  K  K  ,  0^ak<l  (k  =  1,2)  ,  (4) 

where  c  =  -1,  c^=  1,  i  =  fi.  and  are  real  constants  and  u^fx) 
satisfies  the  A/- condition^  near  and  at  c^.  If  u  is  required  to  vanish  at 
x  =  ±  1 ,  the  end  condition 


u(-l)  =  0  and/or  u(l)  =  0  (5) 

is  a  special  case  of  (4)  when  u^.(±l)  =  0  and  u^(x)  satisfies  the  (ji  >  a^)- 
condition.  The  fundamental  function  f  and  the  constraint  functions  f^ 
are  assumed  to  be  at  least  twice  continuously  differentiable  with  respect 
to  their  arguments  u  and  v,  and  continuous  in  x,  but  are  otherwise 


^  A  function  u(x)  is  said  to  satisfy  the  H6lder/^  (ji)-condition  on  path 

L(-l  <  x  <  1)  if,  for  any  two  points  x  ,x  of  L,|u(x)-u  (x  )|^  A|x  -x  |f*  , 

12  2  1  2  1 

where  A  and  [i  are  positive  real  constants,  A  is  called  the  Holder 
constant  and  fX  the  Holder  index. 
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arbitrary.  The  notation  Hx[u]  for  the  finite  Hilbert  transform  of  u,  as 
defined  in  (1),  will  be  used  throughout.  It  may  be  remarked  here  that  the 
solution  of  a  maximum  problem  can  be  deduced  from  this  minimum  one  by 
changing  the  sign  of  the  fundamental  function  in  (2). 

The  Holder  condition on  u(x),  with  specific  index  \i,  and 
the  end  conditions  (4),  or  (5),  are  generally  required  on  physical  grounds. 
Analytically,  it  is  known  (see,  e.g, ,  Muskhelishvili  1953,  §20)  that  if  the 
Holder  index  of  u  is  fx  <  1,  then  v(x)  given  by  (1)  is  also  Holder  con¬ 
tinuous  for  -1  <  x  <  1,  with  fi  <  1.  In  addition,  near  the  end  points 
x  =  ±  1 , 

v(x)=  Tu(?l)log(x±l)  +  v*{x)  (or+ip  =  0,  u(Tl)*0)  ,  (6) 

~  ±  (cot(Q'+iphr)u^(±l)/{xTl)Q'+1P+v^(x)  (oH-ip  £  0)  ,  (7) 

where  v^fx)  satisfies  the  ^/-condition  near  and  at  c  =  -1  or  1  in  the  case 
of  (6)  and  also  in  the  case  of  (7)  if  a  =  0;  however,  for  0  <  a  <  1  in  (7), 

i  ra° 

v^(x)  =  v^^(x)|x  -  c I  ,  with  aQ  <  a  and  v^^(x)  satisfies  the  ^/-condition. 
Thus,  v(±l)  will  be  bounded  either  when  u(*l)  =  0  or  when  a  =  ,  p  =  0. 

If  both  u(±l)  and  v(±l)  are  required  to  be  bounded,  condition  (5)  must  be 
enforced. 

The  original  problem  is  equivalent  to  the  minimization  of  a  new 
functional 

^1 

l[u]  f(u(x),  v(x),x;  .  .  .,  X-M)dx  ,  (8) 

-1 

With  M 

^ (u>  v  >  x>X  .  .  .,  \j^j)  =  fo(u,V,x)“  ^  \  £  [  f|  (u,  V,  x)-C|  ]  ,  (9) 

f  =1 


where  u(x),  v(x)  are  related  by  (1),  X.  ,  .  .  . ,  X  ^  are  undetermined 
Lagrange  multipliers.  We  define  an  admissible  function  as  any  function 
u(x)  which  satisfies  the  Holder  condition  //  (fx  <  1),  the  isopermetric 
constraints  (3),  as  well  as  the  prescribed  end  conditions  (4)  or  (5);  and  we 
define  the  optimal  function  as  an  admissible  function  which  minimizes  the 
functional  l[  u] . 


3. 


The  Necessary  Conditions  of 


Lilli  cULIL> 


Let  u(x)  denote  the  required  optimal  function.  A  function  6u  will 
be  called  an  admissible  variation  if,  for  all  sufficiently  small  positive 
constant  c,  u^(x)  =  u(x)  +  1 6u{x)  is  an  admissible  function.  The  variation 
in  v  which  corresponds  to  an  admissible  variation  6u,  such  that 
v^x)  =  v(x)  +  e6v  =  Hx[u^]  =  Hx[u  +  £^u] .  is  found  from  (1 )  to  be  the  Hilbert 
transform  of  6u, 


6v(x)  =  Hx[6u]  (-l<x<l)  .  (10) 


If  6u  is  an  admissible  variation,  then  l[  u-N  6'^j  is  a  function  of  c  which 
has  an  extreme  value  when  c  =  0. 

The  variation  of  the  functional  I  due  to  the  variations  6u  and 


6v  is 


.1 

I 

AI  =  ^  f(u+«6u,  v+£6v, 

-1 


1 

x)dx  f(u,v,x)dx  , 
-1 


(ID 


For  sufficiently  small  c  ,  expansion  of  the  above  integrand  in  Taylor's 
series  yields 

2  3 

Al  =  1 61  +  -j—  6ZI  +  ■jp  63I  +  .  .  .  ,  (12 ) 

where  the  first  variation  61  and  the  second  variation  6ZI  are 


6l[u,6u]=j  (f^6u  +  fv6v)dx  ,  (13) 

-1 

1 

6Zl[u,6u]  [fuu(6u)2+  2fuv6u6v  +  f^^v)2  ]dx  ,  (14) 

-1 

in  which  the  subindices  denote  partial  differentiations,  and  6v  is  given 
by  (10).  The  variations  61,  6ZI,  .  .  .  depend  on  6u  as  well  as  u.  For 
l[u]  to  be  minimum,  we  must  have  for  all  admissible  variations  6u, 


6l[u,  6u]  =  0 


(15a) 
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and 

62l[u,6u]  ^  0  .  (15b) 

Equation  (15a)  assures  that  I  is  extremal,  and  with  the  inequality  (15b), 

I  is  therefore  a  minimum. 

As  6u  and  6v  are  related  by  (10),  substitution  of  (10)  in  (13) 
reduces  (15a) to 

1 

61  =  j*  (fu-Hx[  fv]  )6u(x)dx  =  0  (16) 

-1 


after  interchanging  the  order  of  integration,  which  is  permissible  (see, 
e.g.,  Tricomi  1957,  §4.  3)  if  the  functions  fv(u,v,x)  and  6u(x)  belong  to 

the  classes  I.  and  L  ,  respectively  (in  the  basic  interval  -1  ^  x^  1), 

Pi  Pz 

and  if 

P,‘+  P2"'«  1  .  (17) 


This  condition  will  be  tacitly  assumed  to  be  satisfied.  Since  6u(x)  is 
otherwise  arbitrary,  the  factor  in  the  parenthesis  of  the  integrand  in  (16) 
must  vanish  identically  for  |x|  1,  giving  the  following  singular  integral 

equation  of  the  Cauchy  type, 

j  f1  f  (u(t),  v(t)J) 

fu(u(x),v(x),x)  =  Hx[fy]  =  Ig:  v-  - -  dt  (|xUl)  .  (18) 

-1 

This  integral  equation  is  analogous  to  the  Euler  differential  equation  in 
the  classical  theory  of  calculus  of  variations  when  the  fundamental  function 
is  of  the  form  f  =  f(y(x),  dy/dx,  x).  Equation  (18)  is  generally  non-linear 
in  u(x)  and  v(x)  unless  f  is  a  polynomial  of  second  degree  in  u  and  v. 
The  extremal  solution  is  determined  by  solving  the  pair  of  coupled  singular 
integral  equations,  (18)  and  (1),  under  conditions  (3)  and  (4).  Satisfaction 
of  (18)  is  a  necessary  condition  for  the  existence  of  an  extremum  of  l[u] . 

We  now  suppose  that  (18),  (1),  (3)  and  (4)  can  be  solved  for  an 
extremal  function  ufxjC^,  .  .  . ,  Cj^),  which  involve  the  constants  of 

constraint  . as  Parameters*  Under  what  condition  does 

this  extremal  solution  satisfy  the  inequality  (15b),  so  that  it  actually  provides 


a  minimum  of  l[u]  ? 


In  order  to  answer  this  question,  we  examine  the  second  variation 

6^.  Consider  the  case  in  which  f  (u,v,x),  f  (u,v,x),  f  (u,v,x)  and 

6u(x)  are  all  Holder  continuous  on  (-1,1),  and  they  belong,  respectively, 

to  the  classes  L  ,  with  p,  so  limited  that  interchange  of  the  order  of 
Pk  k 

the  following  integrations  can  be  justified  (see  (17)  ).  The  second  term  on 
the  right-hand  side  of  (14)  can  then  be  written  as 


.1 

~  ~  2  \  Hx[fuv(u,  v,t)6u(t)]6u(x)dx 
-1 


1  [  £ 

=  v J  T  t  -  x 

-1  -1 


6u(t)6u(x)dtdx 


The  first  step  follows  from  substitution  of  (10),  the  second  step  is  a  result 
of  interchanging  the  order  of  integration,  and  the  last  step  is  the  mean  of 
the  two  preceding  lines.  By  similar  operations,  the  third  term  in  (14) 
can  be  written 


1  1 

j*  fvv(u,v,x)[6v(x)J2dx  =  J  f^u.v.x^tSuttjjH^autsJjdx 
-1  -1 


9 


where,  in  the  last  step,  use  has  been  made  of  the  Poincare- Bertrand 
formula, 


9 


(a  <  x  <  b) 
(19) 


(see,  e.g.,  Muskhelishvili  1953,  §23).  Combining  the  above  results,  (14) 


becomes 

1  1  1 

g(x){fiu(x)}2dx  -  i  j  j  6u(t)6u(x)dtdx 

-1  -1  -1 


(20a) 


where 


g(x)  =  f  (u(x),v(x),x)  +  f  (u(x),v(x),x)  ,  (20b) 

U.U.  V  V 

h(x)  =  fuv(u(x),v(x),x)  +  l^[fvv(u(s)tv(s)f  s)]  .  (20c) 


Since  f  ,f  and  f  are  assumed  to  be  Holder  continuous  in  (-1,1), 
uu*  uv  w 

g(x)  and  h(x)  are  also  Holder  continuous  in  this  interval  (with  the  pos¬ 
sible  exception  of  those  end  points  at  which  f^tu,  v,x)  f  0);  i.e. ,  for  any 


two  points  x  , x  in  the  open  interval  (-1,1), 

1  2 

fi 

lg(x2)  •  gOOl  <  aJx2  -  xj  1  (04  n^<  1,  A^>  0)  ,  (21a) 

U 

|h(x)  -h(x)|  <  A|x  -  x|  2  (0^U<1,  A  >  0)  .  (21b) 

2  1  2  2  1  2  2 

Now  consider  a  special  choice  of  6u(x), 

6u(x)  =  BU(|)  ,  £  =  (x  -  xQ)/«  ,  (22a) 

where  U{£)  is  Holder  continuous  and 

0<U(£K  1  (lei  <  1,  or  |x-xo|<«)  ,  (22b) 

U(|)  =  0  (|el  >  1,  or  |x  -xo|  >€)  .  (22c) 


xq  is  any  fixed  point  in  the  open  interval  (-1,  1),  and  c  is  arbitrarily 
small  so  that  |x  ±  «  |  <  1.  B,  the  bound  of  6u(x),  is  either  positive  or 
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4 


negative,  and  is  chosen  so  small  that  6J  I,  64I,  etc.  ,  can  be  neglected  in 
comparison  with  6ZI.  (For  instance,  B  =  0(€tt),  a  >  0,  will  be  sufficient; 
the  explicit  form  of  U(£)  is  immaterial.)  With  this  choice  of  6u,  (20) 

can  be  written,  by  adding  and  subtracting  a  term,  as 

X  +€  . 

O  1 

=  g(xo)Bz  ]  U2(£)dx  +  R  .  g(xo)Bz«  j*  U2(£)d£  +  R  ,  (23a* 

-1 


x  -t 
o 


where 


X  +€  X  +€ 

R=J  llg(i')-g(xo)]6u(x)  - j  MO  6u(t)dt>6u(x)dx  .  (23b) 

X  -e  X  -€ 

o  o 


Using  the  inequalities  (21),  (22b,  c),  we  obtain  the  bound  for  the  remainder 
R  as 


1 


1 +J1  (i  1  +fi  ^  *  fx  -1 

|rUabz€  lei  ^(ejde+AB2*  zj  j  h-ii  2  u(e)u(q)dedT! 

-i  -i  -i 


2AB2I  1-tti  I  2ABz  \  1-Hi 


(23c) 


In  the  limit  as  «  -*0,  the  first  term  on  the  right-hand  side  of  (23a)  dominates, 
hence  a  necessary  condition  for  6zI  >  0  is  that  g(xQ)  >  0  for  every 
xo€(-l,l),  or 

fuu(u(x),v(x),x)  +fvv(u(x),v(x),x)>  0  (-1<x<1)  .  (24) 


This  condition,  first  derived  by  one  of  the  authors  (AKW  1969),  is 
analogous  to  the  Legendre  condition  in  the  classical  theory  of  the  variational 
calculus. 

Equation  (24  ..s  a  necessary  condition  to  be  satisfied  by  a  minimiz¬ 
ing  function.  A  weaker  form  of  this  necessary  condition  is  f  +  f  ^  0, 

6  1  uu  w  ’ 

meaning  that  the  second  variation  6ZI  >  0  is  still  ensured  if  (f  +  f  ) 
b  uu  w 

varnishes  at  a  discrete  number  of  points  of  the  extremal  function  but  is 
positive  everywhere  else. 

It  may  also  be  noted,  by  ainalogy  with  classical  variational  problems. 
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(see,  e.g.  Courant  and  Hilbert  (1953),  Chap.  IV,  §6),  that  strict  inequality 
in  (24)  is  not  a  sufficient  condition  for  a  minimum.  To  find  a  sufficient 
condition  we  expand  l[u+«6u]  by  Taylor's  theorem  with  a  remainder  after 
two  terms.  Thus, 

2 

I[  u+«6u]  =  l[u]  +  €  6l[  u,  6uj  +  ~y  62l[  u+q€  6u,  6u]  (0  <  q  <  1 ) 

If  u  is  an  extremal  function,  i.  e.  ,  6l[u,6u]  =0,  then 

2 

l[u+€6u]  -  l[u]  =  -j-  6Zl[u+qf6u,  6u]  (0<q<l) 

Now  suppose  that  inequality  (15b)  holds  not  just  for  the  extremal  u,  but 
for  all  admissible  functions  u.  Then  we  may  set  «  =  I  in  the  above 
equation  to  give  the  condition 

l[u+6u]  -  l[u]  =  62l[  u+q6u,  6u]  >0  (0  <  q  <  1)  , 

which  is  sufficient  to  show  that  the  extremal  u  actually  minimizes  I. 

Based  on  the  foregoing  argument,  a  sufficient  condition  for  a 
minimum  is  that  the  quadratic  form,  in  6u  and  6v,  in  the  integral 
representation  (14)  of  62I  be  positive  definite  for  all  admissible  u  and 
6u  (and  hence  all  admissible  v  and  6v  by  (7 )  and  (10)  ),  that  is 


f 

uu 


>  o. 


and 


f  f  -f2  >  0 

uu  vv  uv 


(-1  <  X  <  1) 


(25) 


for  all  admissible  u  and  v.  This  simple  but  rough  sufficient  criterion 
is  a  more  restrictive  inequality  than  (24). 


4.  Quadratic  Functions;  the  Fredholm  Integral  Equation 

The  least  difficult  case  of  the  extremal  problems  in  this  general 
class  is  when  l[u]  is  a  quadratic  functional,  or  when  f  in  (8)  is  quadratic 
in  u  and  v,  since  the  integral  equation  (18)  is  then  linear  in  u  and  v. 

It  is  instructive  to  investigate  this  case  first,  since  the  system  of  singular 
integral  equations  (18)  and  (1)  can  then  be  reduced  to  a  single  Fredholm 
integral  equation  of  the  second  kind,  or,  in  certain  special  cases,  the 
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method  of  singular  integral  equations  can  be  employed  to  obtain  an  analyt¬ 
ical  solution  in  a  closed  form.  These  solutions  provide  a  basis  for 
comparison  with  some  approximate  methods  which  may  have  a  general 
utility  for  more  comj^icated  non-linear  problems.  Moreover,  the  linear 
formulation  can  often  be  used  as  a  first  approximation  of  an  originally  non¬ 
linear  problem  with  appropriate  modification  of  the  isoperimetric  constraints. 

Let  the  functions  f  and  f^  in  (2)  and  (4)  be  given  by 

fQ(u,v,x)  =  aQUZ  +  2bQUv  +  +  2pQU  +  2qQv  ,  (26a) 

fj(u,  v,x)  =  a^u2  +  2b^  uv  +  c^v2  +  2p^u  +  2q^v  (i=l,2,  .  .  .M)  ; 

(26b) 

the  coefficients  a  ,b  ,c  ,  .  .  .  ,  q.  ,  are  known  functions  of  x,  as- 
o’  o  o  ’ 

sumed  to  be  Holder  continuous  (with  index  0  <  p  <  1 )  in  (-1, 1  ).  Then  the 
function  f  in  (8)  becomes 

f(u,v,x)  =  au2  +  2buv  +  cv2  +  2pu  +  2qv  ,  (27) 


where 


M 

\ 


a(x)  =  aQ(x)  -  /  X^ai  (x) 

itl 


etc. 


(28) 


The  integral  equation  (18)  now  reads 


a(x)u(x)  +  b(x)v(x)  +  p(x)  =  Hx[  bu  +  cv  +  q]  ( j  x  |  <  (29) 


The  necessary  condition  (24),  obtained  from  the  consideration  of  the 
second  variation,  becomes 

a(x)  +  c(x)  >  0  ( | x |  <  1 )  ,  (30) 


which  can  be  checked  only  when  the  X 's  in  (28)  are  determined. 

The  coupled  integral  equations  (29)  and  (1),  both  of  the  Cauchy 
type,  can  always  be  reduced,  under  certain  assumptions,  to  a  Fredholm 
integral  equation  of  the  second  kind,  with  a  regular,  symmetric  kernel. 
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The  required  assumptions  are  that  the  coefficients  a,b,  .  .  .  ,q,  as 
well  as  the  solution  u,v,  are  Holder  continuous  on  (-1,1),  and  that  c(x) 
and  u(x)  belong  to  the  classes  L  and  L  ,  respectively,  in  the  closed 

Pi  P2 

interval  [-1,1],  with  p  and  p  satisfying  (17).  In  fact,  substitution  of 
(1 )  in  (29)  yields 

a(x)u(x)  +  b(x)Hx[u]  -  Hx[bu]  -  Hx[  c(t)Ht[u]  ]  =  Hx[q]  -  p(x)  . 


The  second  and  third  terms  on  the  left  side  of  (31)  combine  to  give 

1 

b(x)Hx[u]  -  Hjbu]  =  -  “(t)dt  • 

-1 

Under  the  aforementioned  assumptions,  the  last  term  on  the  left  side  of 
(31)  can  be  rewritten,  using  the  Poincare -Bertrand  formula  (19),  as 

1  1 

-Hx[c(t)Ht[ul  1  =  C(X)U(X)  +  T  )  U(t)dt^  (s--t)(^xr  • 


Thus,  (31)  reduces  to 


{a(x)+c(x)}u(x)  +  j  K(t,x)u(t)dt  =  Hx[q]  -  p(x)  (|x|  <  1) 
-1 


where 


y(t  1L\  .  1  h(t)-b(x)  1  £  c(s)ds 

K(t'x)  -  7  t  -  X  +  7  J  (S-t)(s-J 


(32a) 


(32b) 


If  we  assume  that  condition  (30)  is  satisfied,  we  may  define  a  new  variable 


u(x)  =  (a(x)  +  c(x)}2u(x)  , 


in  terms  of  which  (32a,  b)  become 


~  r  ~  ~  _i 

u(x)  +  j  K(t,x)u(t)dt  =  {a(x)  +  c(x)}  z{Hx[q]  -  p(x)}  (|x|<l)  , 


(34a) 


with 
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K(t,x)  =  {a(t)  +  c(t)}  2  {a(x)  +  c(x)}  2  K(t,  x)  .  (34b) 

This  is  a  Fredholm  integral  equation  of  the  second  kind,  with  a  regular 
symmetric  kernel,  for  which  a  well-developed  theory  is  available.  The 
kernel  and  the  right-hand  side  of  (34a)  will,  in  general,  contain  unknown 
Lagrange  multipliers.  Ideally,  the  integral  equation  can  be  solved  first 
for  arbitrary  values  of  X^'s,  which  can  then  be  determined  by  the  M 
constraints  (4).  Finally,  condition  (30)  should  be  checked. 

5.  Analytical  Solutions  by  the  Method  of  Singular  Integral  Equations 

In  the  general  case  when  the  coefficients  a,  b,  and  c  are 
arbitrary  functions  of  x,  the  solution  of  the  system  of  singular  integral 
equations  (29),  (1)  has  not  been  found  in  closed  form  (for  a  general  discus 
sion,  see  Muskhelishvili  1953,  Part  IV).  How  er,  when  the  coefficients 
a,  b,  c  satisfy  certain  conditions,  the  system  of  equations  (29)  and  (1) 
can  be  reduced  in  succession  to  a  single  singular  integral  equation  of  the 
Carleman  type,  which  can  be  solved  in  turn  by  known  methods,  yielding 
the  final  solution  in  closed  form.  These  analytical  solutions  are  of  great 
interest,  since  in  their  construction  there  are  definite  degrees  of  free¬ 
dom  for  choosing  the  strength  of  the  singularity  of  the  solution  u(x)  at 
the  end  points  x  =±  1.  With  these  possibilities,  the  singular  behavior  of 
u  and  v  near  x  =  ±  1  can  be  explicitly  analysed.  The  following  are 
several  cases  of  fairly  general  interest. 


Case  I 

a(x)  +  c(x)  >  0,  b(x)  =  bQ±  (ac)2,  ( j x|  <  1 )  ,  (35) 

i 

where  bQ  is  a  constant,  and  the  function  (ac)2  stands  for  a  definite 
(say  positive)  branch.  The  first  condition  of  (35)  is  just  the  necessary 
condition  (30)  for  the  optimality,  b(x)  assumes  either  of  the  two  expres¬ 
sions  . 

Multiplying  (1)  by  bQ  and  subtracting  it  from  (29),  we  obtain 


a2<P±(x)  =  Hx[±c2v±]  +4<(x) 


(|x|  <  1) 


(36a) 


15 


where 

®  (x)  =  a2u  +  c2v  ,  4i(x)  =  H  [q]  -  p(x)  .  (36b) 

Z  X 

1 

In  (36),  the  upper  (+)  sign  is  for  b  =  bQ+  (ac)2,  and  the  lower  (-)  sign 
for  b  =  bQ-(ac)2.  We  note  that  the  two  cases  in  (35)  are  the  only  form  of 
b  for  which  (29),  with  the  aid  of  (1),  can  be  reduced  to  an  equation  for  a 
single  variable  <f>+  ot  <j>  .  Now,  (36)  is  a  singular  integral  equation  of 
the  Carleman  type,  the  general  solution  of  which  is  known  (see,  e.g.  , 
Muskhelishvili  1953,  Tricomi  1955).  The  solution  of  (36)  is  found  to  be 


<P±(x)  = 


a24«(x) 

a+c 


c2(tty(t)dt 

Z±(t)(t-x) 


+  pkm^\ 


(|x|  <  1)  (37) 


where  and  are  polynomials  in  x  of  degree  m, 

Z±(x)  =  (1+x)  l(l~x)  2  (a+c)2exp  {±r(x)}  , 

r(x>=2T  Hj^gCtl)]  ,  G(x)  =  (a2  +ic2  )/ (a2  -ic2  )  , 


(38a) 

(38b) 


and  n^,  n^  are  integers  to  be  chosen,  together  with  integer  m,  accord¬ 
ing  to  the  following  rule.  The  function  logG(x)  stands  for  a  definite 
branch  (say  its  principal  branch)  and  is  one-valued  on  [-1,1].  Since  a 
and  c  are  assumed  to  be  Holder -continuous  and  since  a  +  c  >  0  on 
[-1,1],  clearly  G(x)  t  0,  oo.  Consequently,  T(x)  will  be  Holder -con¬ 
tinuous  in  the  open  interval  (-1,1),  but  may  have  logarithmic  singular¬ 
ities  at  the  end  points  x  =  -1,  x  =  1,  unless  logG(x,  )  =  0.  In  general, 

12  X 

r(x)  =  aklog|x  -xkl  +r0(x)  (k  =  1,2)  ,  (39a) 

ak  =  +  ivk  =  (_1  ^  Tin  log  G(xk)  *  (39b) 

where  rQ(x)  remains  bounded  near  and  at  xk_  Hence,  by  (38), 

n  ±<7  n±<7 

Z±(x)  =  w±(x)(l+x)  1  ‘(l-x)  2  2  ,  (40) 

where  w±(x)  are  nonvanishing  functions  in  [-1,1], 
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Now  the  integers  n,  n  are  selected  to  satisfy  either  of  the  con¬ 
ditions 

-l<nk±^ik<0  ,  (41a) 

04  nk±^k<l  (k  =  1,2)  ,  (41b) 

1_ 

in  which  the  upper  (+)  sign  is  for  the  case  b  =  bQ  +  (ac)2  ,  and  the  lower 
(-)  sign  for  b  =  bQ  -  (ac)2.  The  values  of  n,  n^  so  selected  may  be 
different  for  the  two  end  points  (  k  =  1, 2).  Condition  (41a)  is  for  an  end 
point  at  which  Z±(x)  may  admit  a  branch-point  type  (but  integrable) 
singularity,  whereas  condition  (41b)  insures  Z  (x^  =  0,  as  may  be 
required  on  physical  grounds. 

The  selected  n  and  n  determine  the  "end-point  index" 
x  2 

K  =  -(n  +  n  )  ,  (42) 

which  dictates  the  degree  m  of  the  polynomial  f>m(x): 

(i)  For  K  ^  0,  Pm(x)  can  be  taken  an  arbitrary  polynomial  of  degree 

m  4  fC  -  1  ,  if  #C  >  0  ;  and  pm(x)  =  0  »  if  K  -  0  . 

(43) 

(ii)  For  K  <  0,  P^fx)  =  0,  and  the  solution  is  valid  if  and  only  if  the  fol¬ 
lowing  orthogonality  conditions  are  satisfied  by  ^  (x): 

1 

j*  tM(t)[>Mt)/Z±(t)]dt  =  0  (1=0,1,  .  .  .  .-(/c  +  i))  . 

“ 1  (44) 

With  <p±(x)  so  determined,  u(x)  can  be  solved  in  succession  by- 
substituting  (1)  in  (36b),  giving 

a2(x)  u(x)±  c2(x)  Hx[u(t)]  =  <p±(x)  ( | x j  <  1 )  ,  (45) 

which  is  again  a  singular  integral  equation  of  the  Carleman  type.  The 
general  solution  of  (45)  is 
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where  Q^(x)  are  polynomials  in  x  of  degree  m,  Z±  being  the  same  as 

given  by  (38).  In  this  solution,  the  integers  n  ,n  in  the  expression  for 

1  2 

Z  and  the  degree  m  of  are  determined  as  before  Finally,  sub¬ 

stituting  (37)  in  (46),  and  making  use  of  the  Poincare -Bertrand  formula 
(19),  we  obtain  the  solution 


i 

in  which  the  upper  sign  is  for  b  =  bQ  +  (ac)2  ,  and  the  lower  for  b  =  bQ  -  (ac 
The  corresponding  v(x)  is  determined  from  (1)  upon  substitution  of  u(x). 
Note  that  this  solution  remains  valid  even  when  a(x)  and  c(x)  become 
zero  or  Negative  in  [-1,1]  as  long  as  a  +  c  >  0  over  the  entire  interval. 

The  final,  and  perhaps,  most  crucial  step  is  to  examine  the  be¬ 
havior  of  u  and  v  near  the  end  points  x  =  ±  1,  and  to  ascertain  if  certain 
physically  required  end  conditions  can  be  satisfied  by  the  analytical  solu¬ 
tion.  This  task  is  greatly  facilitated  by  the  fact  that  the  solution  is  expres¬ 
sed  in  an  explicit,  closed  form,  as  given  by  (47),  in  this  case.  In  fact, 
taking  the  principal  branch  of  log  G(x),  we  see  that  -ir  <  argG(x)  <  tt 
for  a(x)  +  c(x)  >  0  in  |x|  <  1.  Hence,  by  (39),  we  find  the  bounds 
-  0,  0^M2<  and  by  (40), 

n  ±pi  n  ±pi 

|z±(x)|  =  |u±(x)|(l+x)  1  Ml-x)  2  2  .  (48) 


According  to  rules  (41a,  b),  any  one  of  the  following  three  cases  for  the 
solutions  <p  and  u  are  possible: 

(i)  n  =  0,  n  =  0  //c  =  -(n  +  n  )  =  0).  Then  F*  =  0,  Q*  =  0,  and  Z  Ax)  is 
'  i  2  12  m  ’  m  ’  + 

singular  at  x  =  -1 ,  Z  (x)  singular  at  x  =  1; 


to  I*- 


18 


(ii)  n  =  0,  n  =  -1  for  Z  ,  and  n  =  -1,  n  =  0  for  Z  ,  both  correspond- 

i  2  +  i  *  2  • 

ing  to  K  =  1 thus  permitting  and  to  be  nonvanishing  con¬ 

stants,  with  the  resulting  <p  and  u  both  singular  at  the  two  ends 
x  =  ±  1; 

(iii) n  =1.  n  =0  for  Z,,  and  n  =0,  n  =1  for  Z  ,  both  selections 

12  +  1  -t2± 

corresponding  to  k  -  - 1,  so  that  P*  =  Q*  =  0;  a  solution  which  is 

mm 

regular  at  both  x  =  ±  1  is  possible  if  and  only  if 
1  1 
J  c2(x)»|>(x)/Z±(x)dx  =  0  ,  and  j  <p±(x}/ Z±(x)dx  =  0  .  (49) 

-1  -1 

The  last  condition  can,  of  course,  be  released  if  c(±l)  =  0.  If,  in  addition 
to  (49), 

4*(+i)  =«!»(-!)  =  0  ,  (50) 

then  by  (37)  and  (46),  <p  (±1)  =  0  and  u(±l)  =  0.  The  corresponding 
solution  of  v(x)  will  then  be  bounded  at  x  =  ±  1 . 

Therefore,  existence  of  solutions  of  type  (iii),  with  both  u(x), 
v(x)  Holder -continuous  in  (-1,1)  and  bounded  at  x=±l,  depends  on  the 
satisfaction  of  conditions  (49)  and  (50)  on  a,  c,p  and  q. 

Case  II 


a(x)  >  0  ,  b  =  const.  ,  c  =  0  (51) 

This  is  a  special  limit  of  Case  I;  the  solution  is  easily  deduced  from  (47) 
to  be 

u(x)  =i|»(x)/a(x)  (|x|  <  1)  .  (52) 

This  result  is  also  obvious  from  the  original  equations  (29)  and  (1). 
Boundedness  of  the  solutions  for  u  and  v  now  depends  on  condition  (50) 
only. 

Case  III 


a(x)  =  0 


b  =  const. 


c(x )  >  0  . 


(53) 


This  is  another  special  limit  of  Case  I.  In  this  limit,  logG(x)  =  ii 
hence,  by  (38),  expT  =  (l-x)2(l+x)2.  The  corresponding  solution  for  u 


where 

I  n  +|  n 

X(x)  =  c~2Z  =  (1+x)  1  (1-x)  2  .  (54b) 

This  result  can  be  deduced  from  (47),  or  it  can  be  obtained  directly  from 
(1)  and  (29). 

Case  IV 

a,  b,  c  =  constant  ,  a  +  c  >  0  .  (55) 

This  also  belong  to  Case  I  since  it  is  always  possible  to  find  a 
constant  bQ,  such  that  b  =  bQ  +  (ac)2  when  a,b,  and  c  are  constants. 
Since  G,  as  defined  by  (38b),  is  now  a  constant,  we  have  in  this  case 

i  n  Ta  n  ±o 

Z±(x)  =  (a+c )2(l+x)  1  (1-x)  2  (56a) 

where 

a  =  n  +  iv  =  2^-  log[  (a2  +ic2  )/ (a2  -ic2 )]  .  (56b) 

As  a  is  also  a  constant,  further  simplifications  of  the  expression  for 
the  solution  (46)  and  (47),  are  immediate. 


The  Rayleigh-Ritz  Method 


In  some  physical  problems,  u  and  v  are  required  to  be  Holder  - 
continuous  on  the  closed  interval  [-1,1],  This  means,  in  particular, 
that  condition  (5),  u(tl)  =  0,  must  be  satisfied.  When  the  analytical 
solution  can  be  successfully  obtained,  e.  g.  ,  as  in  Section  5,  the  end  con¬ 
ditions  (5)  can  generally  be  examined  only  at  the  final  stage.  It  is  quite 
possible  that,  on  occasion,  condition  (5)  simply  cannot  be  satisfied,  as 
will  be  illustrated  later  by  examples.  It  might  be  possible  to  remedy 
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this  situation  by  enlarging  the  class  of  acceptable  functions  and  by  modify¬ 
ing  the  isoperimetric  constraints  so  that  a  solution' of  some  physical 
significance  is  reached.  Such  tasks  may  indeed  constitute  the  major 
difficulties  in  the  actual  solution. 

In  case  both  u  and  v  are  required  to  be  Holder  continuous  on  the 
closed  interval  [-1,1],  we  have  seen  that  condition  (5),  u(tl)  =  0,  must 
be  imposed.  An  approximating  procedure  for  this  case  can  be  developed 
by  employing  the  Rayleigh-Ritz  method  (or  the  closely  related  Galerkin 
method),  which  has  been  used  with  great  success  in  other  variational 
problems,  especially  for  numerical  solutions.  This  method  is  applicable 
to  those  variational  problems  which  satisfy  the  sufficient  condition  (24) 
or  the  more  restrictive  condition  (25).  The  central  idea  is  the  construct¬ 
ion  of  a  sequence  of  comparison  functions  that  forms  a  minimizing  sequence. 
A  sequence  u  ,u  ,  .  .  .  ,  Iin,  ...  of  admissible  functions  is  called  a 
minimizing  sequence  for  the  functional  l[u]  if  l[unJ  converges  to  the 
minimum  value  l[u],  where  u  =  limun  as  n  -►  eo  ,  with 

I(u]>4u,]>,.  .  .  £l[u  ]£.  .  .  £l[S]  .  (57) 

12 

whether  l[u]  is  a  minimum  which  is  actually  attained  for  a  function 
u  =  u  remains  an  open  question,  the  answer  to  which  depends,  in  particular, 
on  whether  the  end  conditions  u  kD  =  0  can  eventually  be  satisfied. 

To  construct  a  minimizing  sequence,  we  start  with  a  complete 
system  of  fixed  "coordinate"  functions  <p  (x),  ^(x),  .  .  .  ,  which  are 
Holder  continuous  on  [  -1,1]  and  <?n(±l)  =  0  (n  =  1,2,  .  .  .  ),  so  that 
all  linear  combinations 

un(x>  =  Y^  +  Y  Z<PZ  +  ...  +  V  n<Pn  ,  (58) 

where  y  ,  •  •  .  ,  Yn  are  n  constant  coefficients,  are  admissible 

comparison  functions  for  the  problem.  The  requirement  that  l[un]  =  In 

be  a  minimum  presents  an  ordinary  minimum  problem  for  l[un]  as  a 

function  of  the  n  parameters  y  ,  y  ,  .  .  .  ,  y^,  which  are  determined 

by  the  simultaneous  algebraic  equations  3l[u  ]/3y  =  0  (k  =  1,2,  .  .  .  ,  n), 

n  k 

or  more  explicitly,  when  expressed  in  terms  of  (8), 
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1 

f 


{fu-Hx[fvJ>  ?>k(x)dx  =  0 

-1 


(k=  1,2,  .  .  .  ,n) 


(59) 


Note  that  the  expression  within  the  curly  bracket,  if  set  equal  to  zero  on 

[  -1,1],  yields  just  the  Euler  equation  (18).  This  system  of  equations 

for  Yk'8  can  always  be  fulfilled,  at  least  in  principle  according  to 

Weierstr ass's  theorem,  provided  l[unJ  is  a  continuously  differentiable 

function  of  Yk'8  (we  have  assumed  that  this  is  the  case).  However,  when 

f  and  f  are  nonlinear  or  transcendental  in  u  and  v,  determination 
u  v 

of  the  Yk'8  from  (59),  by  numerical  methods  or  otherwise,  can  be 
extremely  difficult,  particularly  when  n  is  not  small.  Let 
Y  ,  Y,  .  .  .  ,  Yp  be  the  solution  of  (59),  and  set 


Un=I  VkW 


n 

'n 


(60) 


k=l 


k=l 


When  the  fundamental  function  f  satisfies  the  sufficient  condition  (25),  it 
is  obvious  that  l[u  ]  <  l[  u  ].  From  this  it  follows,  by  setting  u  =u 
which  is  always  possible,  that  I[ur]  $  l[un  J,  and  hence  (57)  is  satisfied. 
Therefore,  the  sequence  IT,  u^,  .  .  .  given  by  (60)  is  indeed  a  minimiz¬ 
ing  sequence. 

The  complete  set  of  coordinate  functions  can  be  chosen  in  various 
ways.  Particularly  suitable  for  the  present  class  of  problems  is  the  trigono¬ 
metric  functions  $>n(x)  =  sinnO  (x  =  cos  Q,  n  =  1,2,  .  .  . ),  since  both 

u  and  v  then  become  simply  finite  Fourier  series: 
n  n 

(04  6  =  cos  1  x  4  it)  .  (61) 

This  is  equivalent  to  expanding  u^  and  vn  in  series  of  Tchebichef 
polynomials , 

n  n 

un(x)  =  (1_x2)2  ^Ykuk-i  M  »  vn(x)  =  "XVkTk(x)  (lxl  ^  »  (62a) 

k=l  k=l 


n 

5. -2 


Y^sinkS 


v  = 
n 


k=l 


n 

-l 


Y^cos  k0 


k=l 
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where 

U  (cos  0)  =  sinn0/sin0  ,  T  (cos0)  =  cosn0  (n=l,2,  .  .  .) 
n-i  n 

(62b) 

Tricomi  (19F»)  also  employed  expansions  in  series  of  a  Jacobi  polynomial 
^  )(x)  multiplied  by  its  weighting  function  (l-x^l+x)^  .  Another 
complete  set  of  functions  vanishing  at  x  =  ±1  is  the  polynomials 

=  (l-x2^11  1  (n  =  1,2,  .  .  .);  the  corresponding  expression  for  v  , 
however,  is  more  lengthy. 

For  quadratic  functionals,  the  problem  of  determining  the  Yk's 
reduces  simply  to  a  linear  algebra.  In  fact,  substitution  o£  (27)  and  (58) 
in  (59)  yields 
n 

XKjk7k  =  Lj  (j  =  1,2,  .  .  .  ,  n)  ,  (63a) 

k=l 


where,  with  suitable  rearrangement  of  terms  by  interchanging  the  orders 
of  integration, 

(a^+bfr  .Hj *k]  -^kHx[ *  j ^cHx[ *  j ] HXK]  } ^  ’  (63b) 

-1 


1 

Vi 

-1 


(63c) 


Clearly,  =  K^j-  In  particular,  when  ^(cos  0)  =  sink0,  we  have 

=  ^{asinkO  sinj0-b  sin(k+j)0+c  cosk0  cos  j0}sin0  d0  ,  (64a) 


Lj  =  ^  {q cos  j0-p  sinj0}sin  0  d0 


(64b) 


The  solution  Yk's  of  (63)  contains  M  multipliers  (\^  ,  .  .  .  ,  \^), 
which  can  be  determined,  in  principle  when  n  >  M,  in  terms  of 


C  ,  .  .  .  ,  Cw,  by  the  M  constraints  (3). 
l  M 

Success  of  this  method  in  any  particular  case  will  depend  on  the 
proper  choice  of  coordinate  functions  <p’^  and  on  the  convergence  of 
un-u  as  n -*•<».  Whether  the  minimizing  sequence  itself  converges  to 
the  solution  u  (i.  e.  u  =  u)  is  a  difficult  theoretical  question.  In  many 
cases  this  method  may  still  prove  useful  for  numerical  calculations  even 
though  its  convergence  to  the  *exact  solution  is  unproved.  The  method  will 
be  illustrated  by  examples  and  compared  with  known  analytical  solutions 
in  the  next  section. 


7.  Examples 


The  following  examples  are  selected  to  exhibit  the  main  features 
of  the  optimum  solution.  Several  other  problems  of  physical  significance 
will  be  considered  more  fully  elsewhere. 


Example  1  .  The  functional 


j[u]  =  ^  [-u2(x)  -  2TTx*u(x)v(x)]dx 
-1 

with  u  and  v  related  by  (1),  is  to  be  minimized  subject  to  the  constraint 


(65) 


=  j  u(x)dx  -  1 

-1 


(66) 


Here  we  have  by  comparison  with  (27),  (28), 


a  =  1  ,  b(x)  =  -fix3  ,  c  =  0,  p  =  -  —  \  ,  q  =  0  (67) 


This  set  of  coefficients  does  not  coincide  with  any  of  the  special  cases 
discussed  in  Section  5.  In  this  case,  however,  it  is  still  possible  to 
obtain  a  solution  in  closed  form  since  the  associated  Fredholm  integral 
equation  (32)  is  particularly  simple, 


u(x)  K(t,x)u(t)dt  =  j  \ 


(|x|<l)  , 


(68a) 
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where 

K{t,x)  =  (1?  -x*  )/(t-x)  =  t2  +  tx  4-x2  ,  (68b) 

which  is  a  Pincherle-Goursat  kernel.  Hence  u  -  X.  must  necessarily 
be  a  quadratic  function  of  x.  In  fact,  we  find 

u(x)  =  -  3X2)  ,  X  =  89/45  ,  (69) 

tfis  value  of  X  being  given  by  condition  (66).  The  corresponding  v  is 

v(x)  =  Hx[u]  -  3**)  logjigj-  6*j  .  (70) 

Finally,  upon  substituting  the  solution  (69),  (70)  in  (65),  the  minimum 
value  of  J  is 

JQ  =  J[u]  =  89/90  =  0.9888  .  .  .  .  (71) 

In  this  specific  example  it  is  noteworthy  that  the  logarithmic 
singularities  of  v  at  the  end  points  x  =  ±  1  arise  from  the  fact  that 
u(±l)  £  0.  We  further  note  that  the  values  of  u(±l)  depend  on  the  co¬ 
efficients  listed  in  (67),  whereas  condition  (66)  affects  merely  the  uniform 
scale  of  u(x). 

Although  the  exact  solution  shows  that  v(x)  is  logarithmically 
singular  at  x  =  ±  1,  it  is  nonetheless  significant  to  investigate  how 
accurately  the  minimum  value  Jq  can  be  predicted  by  approximate  solu¬ 
tions  of  u  and  v  that  are  Holder  continuous  on  [-1,1].  Since  u(x)  must 
be  even  in  x,  we  choose  the  n-terms  approximate  solution  as 


n  n 

un(cos  6)  =  y  Yksin(2k-1)0  ,  vn(cos0)  =  -y  ykcos(2k-l )6  . 


(72) 


k=l 


k=l 


Condition  (66)  requires  that  Y  =  2 /tt  for  all  n.  Let  Jn  be  the  minimum 
value  of  J[un]  found  by  minimizing  the  quadratic  functional  J[un]  with 
respect  to  V2»\3»  • 


.y  .  We  obtain  the  following  numerical  results 


shown  in  Table  1.  The  y^'s  in  the  last  row  of  Table  1  are  the  Fourier 
components  of  the  exact  solution 
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y.  =  —  \  u(cos  6)sin(2k-l  )6  dd  , 
k  it  J 

o 

where  u  is  given  by  (69).  These  numerical  results  exhibit  the  increasing 
accuracy  of  Jq  predicted  by  the  Holder  continuous  solutions  un  and  v^, 
as  n  increases.  In  fact,  the  error  of  =  J[un],  by  comparison  with 
the  exact  solution  of  Jq  of  (71),  is  already  as  small  as  0.1%  for  n  =  3. 
Although  the  convergence  of  the  coefficients  y^  to  the  Fourier  components 
of  the  exact  solution  is  less  obvious,  particularly  for  k  large,  plots  of 
u(cos  0)  and  Ujq(cos0),  for  example,  show  very  little  difference  except 
near  the  endpoints  6  =  0,ir(x  =±  1).  The  same  holds  true  for  v(cos  6) 
and  Vjq(cos0),  although  now  the  discrepency  at  the  endpoints  is  more 
pronounced  since  v  in  (70)  is  logarithmically  singular  as  x  -*•*  1. 


Example  2 


The  second  example  is  to  minimize 

1 

J[u]  =  j  {u2(x)  +  kv2(x)}dx 
-1 


subject  to 


u(x)dx  =  1 


(73) 


(74) 


where  k  is  a  real  constant,  k  >  -1,  in  order  that  the  necessary  condition 
(30)  is  satisfied. 

In  this  case,  we  have,  by  comparison  with  (27),  (28),  and  (36b) 


This  belongs  to  case  IV  treated  in  Section  5.  Thus,  by  (56), 

Z±(x)  =  (l+k^l+x)"1"^! -x)1*2*0^  (|x|  <  1)  ,  (76a) 


where 


2' 


a  =  u  =  —  tan  1  k2  (k  >  0)  (76b) 

r  IT 

=  iv=  J-log{[l+(-k)^]/[l-(-kf]}  (-1  <  k  <  0)  .  (76c) 

fair 


A  quick  check  with  the  end  conditions  for  k  >  0  shows  that  conditions 

(49),  (50)  cannot  be  satisfied,  implying  that  the  solution  cannot  be  regular 

at  both  x  =  ,t  1.  In  fact,  by  direct  calculation  from  (47),  with 

P  =  Q  =0,  the  optimum  solution  is  found  as  follows 
m  m 

1  1 

Case  (i)  k  >  0  (tanfiir  =  k2  ,  0  <  fl  <  ^  ) 


u(x)=  cosMu|[^|]  +[|~]  | 


(77a) 

(77b) 


This  solution  is  readily  verified  by  taking  note  of  the  identity 

Hx[(t7f)  /]=±(cotfiir)|T^)  T  csc^  ' 

which  can  be  shown  by  a  contour  integration  of  (t-l/*(t+l)  ^(t-x)  1  en¬ 
circling  the  real  t-axis  from  t  =  -1  to  1  in  the  complex  t-plane  (or  see 
Tricomi  1957,  p.  181).  The  Lagrange  multiplier  X  is  determined  by 
the  constraint  condition  (68)  to  give 


X  =  (tan  [in )/  fiir  , 

and  the  corresponding  minimum  value  of  J  is 

JQ  =  \  X  =  (tanfiir)/ {Zfiir)  (0  <  p  <  ~  ) 

i 

Case  (ii)  -1  <  k  <  0  (tanh  vrr  =  ( -k)2  ,  v  >  0)  . 

u(x)  =  2"  X  cosh  VTT  COS  |  V  log  J  , 
v(x)  =  2"  X  cosh  vtt  coth  vir  sin |  v  log  j  , 


(78a) 


(78b) 


(79a) 

(79b) 


with 
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X  =  (tanh  vtt)/  vtt 


J  =  *-  X  =  (tanh  vtt)/ (2  vtt) 

O  b 


(80a) 


(80b) 


The  above  solution  can  be  deduced  from  Case  (i)  by  analytic  continuation 
of  the  parameter  [i  =  iv,  or  it  can  be  verified  by  making  use  of  the  con¬ 
version  formulae 

Hx|^sin|v  log  =  “  coth  vtt  cos  |  v  log  y^-|+  csch  vn  , 

Hx  jcosjv  log  y^-)J=  coth  vtt  sin  Jv  log  , 

which  can  be  shown  by  a  contour  integration  of  (t-l)lv(t+l)  lv(t-x)  1 
circumventing  the  real  axis  from  t  =  -1  to  1  in  the  complex  t-plane. 

A  third  limiting  case  is  k  =  0,  or  c  =  0,  in  which  case  it  fol¬ 
lows  from  (52)  that  u(x)  =  i-  ,  X  =  1  and  Jq  =  j  .  This  result  agrees 

with  Cases  (i)  and  (ii)  in  the  limit  as  k  0. 

To  compare  this  solution  with  the  Rayleigh-Ritz  method,  we  choose 
k  =  -  1/4,  in  which  case 

Jq  =  1/(2  log3)~  0.455119  . 

The  discretized  Fourier  method,  again  using  the  expression  (72)^ gives 
the  results  shown  in  Table  2.  The  apparent  convergence  of  the  Fourier 
series  method  to  the  exact  solution  is  again  exhibited. 
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THEORY  OF  OPTIMUM  SHAPES  IN  FREE -SURFACE  FLOWS 

PART  n. 

Optimum  Profile  of  Spray  less  Planing  Surface 

by 

T.  Yao-tsu  Wu,  Arthur  K.  Whitney 

California  Institute  of  Technology 
Pasadena,  California 

The  purpose  of  this  work  is  to  evaluate  the  optimum  profile  of  a  two- 
dimensional  plate  producing  the  maximum  hydrodynamic  lift  while  planing 
on  a  water  surface,  under  the  condition  of  no  spray  formation  and  no 
gravitational  effect,  the  latter  assumption  serving  as  a  good  approximation 
for  operations  at  large  Froude  numbers.  By  employing  a  recently  develop¬ 
ed  theory  of  variational  calculus  involving  singular  integral  equations,  the 
lift  of  die  sprayless  planing  surface  is  maximized  under  the  isoperimetric 
constraints  of  fixed  chord  length  and  fixed  wetted  arc -length  of  the  plate. 
Consideration  of  the  extremization  yields,  as  the  Euler  equation,  a  pair 
of  coupled  nonlinear  singular  integral  equations  of  the  Cauchy  type,  which 
are  subsequently  linearized  to  facilitate  further  analysis.  The  analytical 
solution  of  the  linearized  problem  has  a  branch-type  singularity,  in  both 
pressure  and  flow  angle,  at  the  two  ends  of  plate.  In  a  special  limit, 
this  singularity  changes  its  type,  emerging  into  a  logarithmic  one,  which 
is  the  weakest  type  possible.  Guided  by  this  analytic  solution  of  the 
linearized  theory,  approximate  solutions  have  been  calculated  for  the 
nonlinear  problem  using  the  method  of  discretized  Fourier  expansions, 
and  the  numerical  results  compared  with  the  linearized  theory. 
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PART  II 

Optimum  Profile  of  Sprayless  Planing  Surface 
1.  Introduction 

The  problem  of  planing  surface  has  received  much  interest  in  the 
past  as  a  device  for  producing  hydrodynamic  lift,  while  moving  forward  on 
a  water  surface.  Most  of  the  early  theoretical  studies  were  based  on  the 
linearized  theory,  taking  into  account  the  effect  of  gravity  for  the  range  of 
moderate  to  large  Froude  numbers,  and  assuming  that  the  spray  sheet  at 
the  leading  edge  of  the  plate  is  thrown  backward  in  the  upstream  direction. 
The  hydrodynamic  drag  on  the  planing  surface  therefore  consists  essential¬ 
ly  of  two  components,  one  due  to  spray  formation  and  the  other  due  to  wave- 
making,  aside  from  the  viscous  skin -frictional  drag,  which  is  generally 
small.  An  exhaustive  survey  of  the  literature  on  the  linear  theory  of 
planing  surfaces  has  been  given  by  Wehausen  and  Laitone  (I960).  A 
crucial  limitation  of  the  linear  theory,  which  seems  to  have  escaped 
proper  recognition,  is  that  the  plate  draft  (or  the  height  of  the  plate 
above,  or  below,  the  undisturbed  water  surface)  cannot  be  arbitrarily 
prescribed.  Loss  of  this  degree  of  freedom  may  be  attributed  to  the 
preassigned  direction  of  the  spray  sheet.  This  limitation  was  removed 
by  Rispin  (1967)  and  Wu  (1967),  who  developed  a  nonlinear  theory  based  on 
the  singular  perturbation  method. 

Of  all  the  previous  investigations,  an  important  contribution  by 
Cumber  batch  (1958)  may  be  singled  out  as  the  only  case  in  which  the  pos¬ 
sibility  was  explored  for  a  planing  surface  to  operate,  at  a  given  Froude 
number,  without  spray  formation  -  -  the  so-called  "smooth  entry"  con¬ 
dition.  This  state  of  operation  immediately  opens  up  the  possibility  of 
further  drag  reduction  by  eliminating  the  spray,  thereby  improving  the 
hydromechanical  efficiency  of  planing  surface. 

This  paper  seeks  to  determine  the  optimum  profile  of  a  two- 
dimensional  plate,  moving  along  the  free  surface  of  an  otherwise  un¬ 
disturbed  water,  without  forming  a  spray  sheet  at  the  leading  edge,  such 
that  for  given  chord  length  and  wetted  arc -length  of  the  plate,  this  profile 
will  maximize  the  lift.  For  simplicity,  the  Froude  number  is  assumed  to 
be  so  large  that  the  gravity  effect  may  be  neglected  as  the  first  approximation,  or 
can  be  evaluated  separately  in  ahigher  order  theory.  The  flow  is  further  assumed 


to  be  inviscid  and  irrotational.  Consequently,  in  the  absence  of  the 
gravitational  and  viscous  effects  as  well  as  spray  formation,  the  planing 
surface  will  encounter  no  drag,  leaving  the  lift  as  the  only  component  of 
the  hydrodynamic  force. 

Aside  for  its  practical  value  in  engineering  applications,  this 
problem  was  selected  originally  as  one  of  the  simplest  in  the  general  theory 
of  optimum  shapes  involving  free  surface  flows,  a  theory  which  may 
have  a  far-reaching  significance  in  its  mathematical  context.  Generally 
speaking,  in  this  class  of  variational  problems,  the  functional  subject  to 
extremization  contains  unknown  argument  functions  which  are  related  to 
each  other  by  integral  equations  -  -  a  consequence  of  the  very  nature  of 
the  mixed-type  boundary  problems.  In  particular,  for  two-dimensional 
potential  flows,  the  integral  equation  is  singular,  of  the  Cauchy  type. 

Thus,  this  situation  is  in  sharp  contrast  to  classical  variational  calculus, 
in  which  the  unknown  argument  functions  are  related  by  differential  equa¬ 
tions.  Consequently,  the  Euler  equation  which  results  from  the  considera¬ 
tion  of  extremization  turns  out  to  be,  in  general,  a  nonlinear,  singular 
integral  equation.  Since  the  methods  of  solution  of  this  equation  are  very 
limited,  more  powerful  methods  are  very  much  desired.  A  preliminary 
mathematical  study  of  this  new  class  of  variational  problems  has  been 
carried  out  by  the  authors  (see  Part  I).  Following  the  same  approach, 
the  present  problem  will  be  investigated  to  provide  useful  solution  of 
hydromechanical  interest.  It  is  hoped  that  this  study  will  stimulate 
further  interest  in  the  development  of  the  general  theory,  and,in  turn,  aid 
the  resolution  of  numerous  fluid  mechanical  problems  of  potential  useful¬ 
ness. 


In  order  to  prepare  for  the  formulation  of  the  optimum  shape 
problem,  we  begin  with  a  consideration  of  the  entire  class  of  two-dimen¬ 
sional  plates,  planing  on  a  water  surface,  which  is  otherwise  undisturbed, 
with  the  plate  profiles  so  adjusted  that  the  entry  of  water  at  the  leading 
edge  of  plate  is  "smooth,"  i.e.  ,  without  forming  a  spray  sheet,  as  shown 
in  Fig.  1.  It  is  convenient  to  choose  the  body  frame  of  reference  so  that 
the  free  stream  velocity  is  U,  in  the  positive  x-direction.  The  resulting 
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flow  is  assumed  to  be  incompressible  and  irrctational.  TheFroude  number, 
Fr  =  U/{g£  )2 ,  based  on  the  chord  length  f  and  the  gravitational  constant 
g,  is  taken  to  be  sufficiently  large  so  that  the  effect  of  gravity  may  be 
neglected.  This  class  of  flows  is  thus  characterized  by  having  no  spray 
sheet  and  no  stagnation  point  inside  and  on  the  flow  boundary,  provided 
the  plate  has  a  continuous  slope. 

By  a  suitable  choice  of  the  origin  and  magnification,  the  complex 
potential  f  =  <p  +  i4»,  <p  being  the  velocity  potential  and  ^  the  stream 
function  ,  is  mapped  onto  the  lower  half  of  a  parametric  l,  =  £  +  iq  plane 
by 

f  =  AU4  ,  (1) 

where  A,  a  real  positive  constant,  is  chosen  so  that  the  plate  is  mapped 
onto  11=  0,  -  1  <  £  <  1,  and  die  free  surface  onto  ti  =  0,  |£  |  >  1.  The 
physical  plane  will  be  denoted  by  z  =  x  +  iy,  in  which  the  x  and  y 
components  of  the  flow  velocity  are  u  and  v,  respectively.  In  terms  of 
the  complex  velocity 


w  =  df/dz  =  u 


.  -i0 

iv  =  qe 


q  =  (u2  +  v2)* 


or  in  terms  of  the  logarithmic  hodograph  variable 


e 


=  tain  (v/u) 


9 


(2) 


w  =  log(U/w)  =  r  +  id  ,  T  -  log(U/q)  ,  (3) 

the  Bernoulli  equation  reads 

P  -  PQ  =  j  P(U2  -  q2)  =  j  pU2(l  -e"Z  T)  ,  (4) 

p  being  the  pressure,  pQ  its  free  stream  value,  and  p  the  fluid  density. 
On  the  free  surface  =  0),  p  =  pQ,  hence 

t'(£)  =  t{£,0-)  =  0  (lei  >  1)  .  (5) 

On  the  plate,  we  denote  the  boundary  value  of  <»>(£,)  by 

u"(£-i0)  =  T_(£)  +  ie"(£)  =  F(S)  +  ip(£)  (lei  <  1)  .  (6) 
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As  for  the  boundary  condition  on  the  plate,  the  simplest  approach  is  to  con¬ 
sider  the  "inverse  problem"  by  prescribing  either  T(£),  or  (3(1),  as  a 
known  function  of  £,  together  with  certain  conditions  to  be  specified  be- 
low.  When  T(£)  is  prescribed,  it  is  required  to  be  Holder -continous  , 
non-negative(to  insure  that  the  pressure  on  the  plate  is  nowhere  less  than 
Pq),  and  T  (£ )  is  required  to  be  continuous  across  the  two  ends  of  the 
plate  at  %  -  ±  1 ,  that  is 

T-(e)  =  r(e)^o  deUn  ,  m 

r(i)  =  r(-i)  =  o  .  (8) 

Under  the  present  assumptions,  the  free  stream  condition  is  simply 

u-0  as  IU  -oo  ,  0  .  (9) 

The  solution  of  the  Dirichlet  problem  prescribed  by  (5),  (7)  and  (9)  is 

1 

«(U=^j  h=ImU0)  .  (10) 

-1 

As  f,  approaches  an  arbitrary  point  t,  -  £  -  iO  on  the  plate,  use  of  the 
Plemelj  formula  (cf.  e.g.  Muskhelishvili  (1953)  )  shows  that  the  real 
part  of  (10)  reduces  to  an  identity,  and  its  imaginary  part  gives 

1 

P(S)  =  TT^  5 *Vr{t)]  <le I  <  1)  ,  (11) 

-i 

in  which  the  symbol  C  over  the  integral  signifies  the  Cauchy  principal 
part  of  the  integral,  and  the  symbol  H|[r]  denotes  the  finite  Hilbert 
transform  of  T  on  [  -1,1]  .  It  is  noted  that  if  T( |)  is  Holder-continuous 
on  [-1,1]  and  if  the  end  conditions  (8)  are  satisfied,  then  P(£),  given 
by  (11),  is  also  Holder -continuous  on  [  -1,1]  (see  Muskhelishvili  (1953) 
§19,29). 

If,  on  the  other  hand,  0  (£)  =  (3(£),  instead  of  t  (£),  is  prescribed 

r($)  is  said  to  be  Holder -continuous  on  [-1,1]  if  for  any  two  points 

!»£  on  [-1,1],  |r(£  )-r(|  )|  <  b||  with  the  Holder  constant  B  >  0, 

1  2  12  12 

and  the  Holder  index  [i  satisfying  0  <  fi  £  1 . 
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for  1 1 1  <  1 ,  one  may  either  solve  this  Riemann -Hilbert  problem  directly, 
or  regard  (11)  as  an  integral  equation  for  T(£).  The  solution  for  T(£) 
satisfying  conditions  (8)  and  (9)  is  found  to  be 

J 

r(£)=  -  I  (l-|2f  £  - £  (;)dt-  ■  (lei  <  1)  .  (12) 

-1  (1-t2)2  (t-|) 

provided  P(|)  further  satisfies  the  orthogonality  condition 

1 

j  P(i)(i-e2)"2de  =  o  .  d3) 

-i 

i 

In  the  above  and  henceforth,  the  function  (42  -l)2  is  defined  to  be  one¬ 
valued  in  the  entire  complex  £  -plane,  cut  from  4  =  -1  to  1  along  the  real 
4 -axis,  so  that  (^2-l)2-,,4  as  |4|  "*■»  for  all  arg4.  Thus, 

(42 -1  j2  “*■*  i(l -42)2  as  4  ±  iO,  |4|  <  1.  It  may  be  remarked  here 

that  the  3ystem  (8)  and  (11)  is  equivalent  to  the  system  (12)  and  (13),  since 
p(|)  given  by  (11)  with  T(|)  subject  to  condition  (8)  satisfies  (12)  and 
(13),  and  conversely,  r(|)  given  by  (12)  with  P(|)  subject  to  (13)  satis¬ 
fies  (11)  and  (8).  Furthermore,  by  virtue  of  condition  (5),  <*>(4)  can  be 
continued  analytically  into  the  upper  half  4  -plane  by 

«(D  =  -  .  (i4) 

The  physical  plane  is  obtained  by  integration  of  w  =  df/dz, 

z(4)  =  AUj  d4/w(4)  =  A  \  eW(;,d4  ,  (15) 

-1  -1 

z(-l)  being  chosen  to  be  the  origin.  The  chord  of  the  plate,  f  ,  and  its 
angle  of  incidence  to  the  free  stream  ^(positive  in  the  clockwise  sense) 
are  given  by 

^1 

I  =  A  j  er(e,cos[p(|)+a]d|  ,  (16) 

-1 

1 

y  er^sin[  p  (|)+o]d|  =  0 
-1 


(17) 
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Equation  (16)  determines  the  factor  A  in  the  transformation  (1),  and 
(17)  states  that  the  angle  of  attack  a  is  referred  to  the  chord  of  the  plate. 
Finally,  the  total  arc -length  S  of  the  plate  is 

1 

S  =  Aj  er^>de  .  (18) 

-1 


The  total  force  F  =  D  +  iL,  D  being  the  drag  and  L  the  lift  act¬ 
ing  on  the  plate,  can  be  determined  from 


B  1 

n  *. 

F  =  i  j  (P-P0Hz  =  j  ip  j  (U2  -w)  -  ^  dC, 

A  -1 

=  j  ipU2A  j  £eW(^-e~W^  j  d£  =  j  ipU2A(^)  eW(^)d£(> 
-1 


(19) 


where  the  contour  of  the  last  integral  encircles  the  plate  counter-clockwise 
in  the  £  -plane  upon  using  the  analytic  continuation  (14).  By  expanding 
«(£),  given  by  (10),  for  large  |^|  ,  we  obtain,  by  the  theorem  of  residues, 

1 

L  =  pU2A\  r(£) d£  ,  D  =  0  (20) 

-1 


Thus,  the  drag  D  =  0,  as  should  be  expected  since  there  is  no  mechanism 
for  producing  drag,  wavemaking  or  otherwise,  in  this  idealized  case. 


3.  The  Optimum  Shape  Problem 

We  now  consider  the  optimum  shape  problem:  In  the  class  of 
functions  T(£)  which  are  Holder -continuous  on  [-1,1],  satisfy  the  in¬ 
equality  condition  (7),  and  the  homogeneous  end  conditions  (8),  find  the 
extremal  arcs  Fq(£)  and  its  conjugate  PQ(£),  mutually  related  by  (1 1 ), 
which  maximize  the  lift  L  under  the  isoperimetric  constraints  of  fixed 
chord  f  and  total  arc-length  S. 

In  what  follows  we  shall  assume  that  the  extremal  arc  has  the 
property  a  =  0  (zero  incidence  of  the  chord)  and  the  symmetry 
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n-S)  =  r(£)  ,  p(-£)  =  -P(S)  •  (21) 

The  fact  that  the  solution,  if  unique,  must  have  this  property  may  be  seen 
by  observing  that  the  extremal  arc  will  remain  extremal  when  the  flow 
direction  is  reversed.  (This  statement  can  actually  be  proved  mathematical¬ 
ly  by  using  a  reverse  flow  argument  and  requiring  that  both  T  and  p  be 
bounded  at  the  two  ends  of  the  plate.  )  Under  this  condition,  (17)  is  then 
automatically  satisfied. 

The  problem  of  maximizing  the  lift  L  =  pU2  L  (see  (20)  )  under 
the  isoperimetric  constraints  of  fixed  chord  i  (see  (16),  now  with  or  =  0) 
and  given  arc-length  S  (see  (18)  )  is  equivalent  to  that  of  finding  the  pair 
of  extremal  arcs  r(£ ),  p  (£  ),  which,  while  satisfying  (7),  (8)  and  (11), 
will  also  minimize  the  new  functional 

1 

l[l\p;A]  =  XI  +  XS  -  L*  =  A  \  f(r(£),p(£);  X  ,\  )d£  ,  (22a) 

12  J  12 

-I 

with  the  fundamental  function  given  by 

f(r,p;X  ,X  )  =  X  er^,cosp(g)  +  X  er(^-  T(£)  (|$|<1)  .  (22b) 


Here,  X  , X 
tive  sign  to 


are  undetermined  multipliers,  and  we  have  assigned  a  nega- 
L  so  that  minimization  of  l[  T,  p  ;A]  corresponds  to 


maximization  of  L  ,  It  is  necessary  to  include  the  coefficient  A  in  the 


arguments  of  I,  since  for  fixed  f  and  S,  A  is  a  functional  of  r  and 


P. 


The  general  variational  problem  of  this  kind  has  been  discussed 
recently  by  the  authors  (see  Part  I).  For  the  problem  at  hand,  the  method 
of  solution  will  follow  the  same  approach,  though  some  modifications  are 
required.  Let  the  set  [  r(| ),  p  (|  );A]  denote  the  optimal  solution  and  let 
[r(£),p(£);A]  be  an  arbitrary  neighboring  admissible  set,  which, by 
definition,  satisfies  (11),  conditions  (7)  and  (8),  and  the  Holder-continuity 
condition.  The  differences  6r  =  r{£)  -  T(£),  6p  =  p  (|)  -  p(£), 

6 A  =  A^-A  form  a  set  of  arbitrarily  small  variations.  Suppose  6r(£)  is 
taken  to  be  a  small,  arbitrary  function  of  then,  since  both  [  F,  p  ]  and 
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[T,p  ]  satisfy  (11),  6p(£)  is  given  by  the  Hilbert  transform  of  6 T, 

6p(e)  =  H£[6r]  (lei  <  1)  •  (23) 

s 


The  variation  6 A,  however,  is  arbitrary. 

The  variation  of  the  functional  I  due  to  the  variations 
[  61",  6(3 , 6  A]  is 

J  1 

AI  =  (A+6  A)  ^  f(r+6r,  p  +6p  )d£  -  A  £  f(r,  p  )d£  . 


-1 


-1 


Expansion  of  the  above  expression  for  sufficiently  small  |6r|  ,  |6p  |  ,  |6A| 
yields 

AI  =  61  +  j  62I  +  63I  +  .  .  .  , 


where  the  first  variation  61  and  the  second  variation  6ZI  are 


61  =  (6  A)  j  f(r,p  )d£  +  A  y  (fr6r  +  fp &p  )d£  , 

6ZI  =  2(6A)j'(fr6r+fp6p)de+Aj[frr46r)2+2frp6r6p+fpp(6p)2]d|  , 


in  which  the  subindices  denote  partial  differentiations,  and  all  integrals 
are  from  £  =  -1  to  1.  For  l[  T,  p  ;A]  to  be  minimum,  we  must  have 
61  =  0  and  62I  >  0  for  arbitrary  6r  and  6A.  From  61  =  0  it  then  fol¬ 
lows  that  the  two  integrals  in  the  expression  for  61  must  vanish  separately. 
The  first  integral  vanishes  if,  by  (22a), 

L^rX^+X^  ,  (24a) 

or  explicitly, 

)  r<£)d£=j  er(e)[\iCosp(£)  +  \2]d£  ,  (24b) 

-1  -1 

which  provides  one  condition  for  the  constant  multipliers  X  ,  X^.  (Note 
that  the  positive  coefficient  A  drops  out  in  (24b).)  For  the  second  integral, 
we  substitute  (23)  for  6(3  ,  then  interchange  the  order  of  integration,  which 
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is  permissible  under  certain  conditions  {see,  e.g.,  Tricomi  1957,  §4.3), 
giving 

i  <fr«r  +  f  6(J)d£  =]  (fp-HgtfpDfingJdg  =  0  .  (25) 

-1  -1 

Since  6F(£)  is  arbitr  ary,  we  obtain  the  following  nonlinear  singular 
integral  equation  of  the  Cauchy  type: 

i  C  fR<r(t),p(t) ) 

fr(r(S),P(S) )  =  H^tfp]  =  i  J  -  dt  (U I  <  i)  , 

-l 

(26a) 

where,  by  (22b) 


^  =  er^>[Xcosp(e)  +  X2]-l  ,  =  -X  er^sinp(£)  . 

(26b) 

For  the  extremal  solution,  (26)  is  to  be  solved  together  with  (11),  as  a 

pair  of  singular  integral  equations  for  F(£)  and  p(£),  subject  to  the 

homogeneous  end  conditions  (8)  and  the  inequality  condition  (7).  The 

extremal  solution,  r(£;X  ,X  )  and  P(|;X  ,X  ),  when  determined  in  this 

12  12 

manner,  will  involve  the  two  constant  multipliers  X  and  X  ,  which  can 

1  2 

be  determined,  most  conveniently,  by  applying  condition  (24),  and  by 
giving  a  specified  ratio  of  the  arc-length  S  to  the  chord  f  ,  say 


S/i  =  1  +K  (K  >  0)  ,  (27a) 


or, by  using  (16),  (18), 

1  1 

|  er(e)d£  =  (1+k  )  J  er(^)cosp(e)de  .  (27b) 

-1  -1 

Since  the  coefficient  A  does  not  appear  in  either  of  the  isoperimetric 
conditions  (-4)  and  (27),  the  problem  of  determining  the  unknown  A  is 
curtailed  altogether.  Finally,  the  optimum  lift  coefficient,  upon  using 
(24a)  and  (27a),  can  be  expressed  as 

CL  =  l/(J.  pU2i)  =  2L*/i  =  2Xt  +  2X2(14K  )  .  (28) 
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This  optimum  lift  coefficient  will  be  a  maximum  if  the  second 
variation  of  I  satisfies  the  inequality  condition  62I  >  0,  which  is  reduced 
to 

1 

-» 

\  [frr46r)2+  2frp6r6p  +f(3p(6p)2]dS  >  0 

-1 

upon  incorporating  (25)  and  noting  that  A  >  0.  A  necessary  condition  for 
the  above  inequality  to  hold  has  been  found  by  Whitney  (1969)  to  be 

frr(r(e),p(e)Hfpp(r(§),p(e)  )>  o  (-1  < 6  <  i)  .  (29) 

The  procedure  for  obtaining  this  result  is  to  first  substitute  (23)forthe  6p 's 
in  the  integrand,  interchange  the  order  of  integration  by  the  Poincare- 
Bertrand  formula,  and,  finally,  to  consider  a  special  choice  of  5T(£) 
which  vanishes  everywhere  except  on  an  infinitesimal  stretch  in  (-1  <£  <  1). 
For  the  present  problem,  with  f  given  by  (22b),  (29)  gives 

frr +  fp p  =  ^2^^^  ®  ’  or  simP1V  X^>  0  »  (30) 

since  the  optimum  solution  P  is  real.  It  may  be  remarked  here  that 
in  (28)  will  be  a  minimum  when  X  <  0.  Condition  (30)  also  shows  the 
importance  of  including  the  arc-length  S  as  a  constraint;  otherwise,  the 
necessary  condition  (29)  cannot  be  satisfied,  and  the  consideration  of 
optimality  must  necessarily  proceed  to  higher  order  varia.ti.ons  of  the 
functional  I,  to  say  the  least. 

The  exact  solution  of  this  problem  is  exceedingly  difficult  for 
several  reasons.  First  of  all,  (11)  and  (26)  are  a  system  of  nonlinear 
singular  integral  equations,  with  a  Cauchy  kernel,  which  have  no  known 
general  method  of  solution.  Second,  it  appears  to  be  very  difficult  to 
incorporate  automatically  the  inequality  condition  (7),  T(^)^  0  for 
| £  |  <  1,  into  the  analysis,  the  only  alternative  being  to  verify  its  validity 
if  and  when  all  possible  solutions  for  T  have  been  obtained.  Furthermore, 
there  is  no  assurance  that  the  homogeneous  end  conditions  (8), 
r(l)  =  T(-l)  =  0,  can  always  be  satisfied.  Finally,  even  when  the  solution 
of  r  satisfying  all  these  conditions  can  be  obtained,  the  determination  of 
the  multipliers  from  (24),  (27)  will  involve  equations  which  are 
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highly  transcendental.  The  foregoing  observations  should  indicate  that  any 
plausible  method  of  solution  by  numerical  iterations  would  most  likely 
meet  great  resistance. 


However,  important  information  about  the  solution  can  be  obtained 
from  the  corresponding  linearized  theory,  which  we  proceed  to  consider 
in  the  following. 


4.  The  Linearized  Theory  (for  (S-i  )/i  «  1)  • 

The  linearized  theory  is  expected  to  provide  a  vSi  first  order 
solution  to  P(£)  and  P(£)  when  the  arc  length  S  is  oiuy  slightly  greater 
than  the  chord  f  ,  or 


S/i  =  1  +  K  ,  0  <  K  «  1  .  (31) 


In  this  limit,  T(£ )  and  p(£)  are  anticipated  to  be  almost  everywhere 
small  on  (-1  <  £  <  1 ),  except  possibly  near  the  end  points  £  =  ±1.  Thus, 
upon  expanding  fp  and  fp  for  small  |  rj  and  |  p  |  ,  and  keeping  only  the 
linear  terms,  Eq.  (26)  reduces  to 

aT(£ )  =  c  H^[  p  (t )]  +  (1  -a)  (|||<1)  ,  (32a) 

where 

a  =  Xt+  \2  ,  c  =  -\i  .  (32b) 

The  linear  system  of  singular  integral  equations  (32)  and  (11)  belongs  to 
the  class  investigated  previously  by  the  authors  (see  Part  I,  Sect.  5),  a 
class  which  cam  be  uncoupled  to  yield  a  set  of  singular  integral  equations 
of  the  Carleman  type  and  then  solved  by  known  methods.  Without  going 
through  the  detailed  analysis,  we  give  below  the  final  solution,  which  can 
be  readily  verified.  The  solution  has  two  branches  according  as  the  co¬ 
efficient 

cr  =  c/a  =  -X  /(X  +  X  )  >  0  or  <0  .  (33) 

1  I  2 

Case  (i)  a  >  0.  The  ranges  of  X  and  X  in  this  case  are  either 
-  1  2 

X  >-X  >  0  or  X  <  -  X  <  0  ;  (34) 

2  1  2  1 

the  first,  according  to  the  necessary  condition  (30),  corresponds  to  the 
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maximum  lift,  whereas  the  second,  to  the  minimum  lift.  In  terms  of  the 
new  parameters,  y  and  fi,  defined  by 

i 

Y=(l-a)/a  and  tanftir  =  &  2  (0  <  fi  <  1/2)  , 

(35) 


the  solution  is  given  by 

r{|)  =  ^  cosfiir  +  (t^T)  j  ' 

p(g)  =  J  cosfiir  cot(Hrp^j||  -  (^j  j  . 


(36a) 

(36b) 


This  solution  is  readily  verified  by  making  use  of  the  formula 


-.[Nil  =  ±  cos  flTT  +  (cotftir) 

which  can  be  derived  directly  by  contour  integration  of  (t-l/^t+l)  ^(t-£)  1 
in  the  complex  t-plane  (or  see  Tricomi  (1957),  p.  181).  We  note  that  the 
above  solution  satisfies  the  inequality  condition  (7),  i.e.  T(£)>  0  dei  <  i), 
but  is  singular,  with  a  branch-type  singularity,  at  the  end  points  |  =  ±1, 
and  thus  fails  to  satisfy  the  homogeneous  end  condition  (8).  This  singular 
feature  of  the  solution  is  perhaps  reasonable,  since  linearization  general¬ 
ly  introduces  singularities  at  those  points  where  the  assumptions  of 
linearization  are  violated.  Whether  these  singularities  can  be  removed 
by  including  the  nonlinear  terms  remains  to  be  seen. 


del  <  i) 


We  now  determine  the  Lagrange  multipliers  \  from  (24)  and 

(27).  To  be  consistent  within  the  framework  of  the  linearized  theory,  all 
the  nonlinear  functions  in  the  integrands  of  (24b)  and  (27b)  will  be  expanded 
for  small  |r|  and  |p  |  up  to  the  quadratic  terms,  because  the  linear 
integral  equation  (32)  actually  follows  from  (26)  by  expanding  the  funda¬ 
mental  function  f(T,  p  )  up  to  the  terms  with  T2,  Tp,  and  p2  .  Thus, 
(24b)  c  d  (27b)  reduce,  sifter  some  rearrangement,  respectively  to 


1  1  1 

2-yj*  r(e)d£  +  ^y  r2(£)d£  +  Ia^  p2(£)d£=0  ,  (37) 

"-1  -1  -1 

r1  r1  1 

2  t  j ^  r (£ )d£  +  \  J  r2(£)d£  -  !+*-  J  p2(£)d£  =  0  , 


44 


45 


Substituting  these  expressions  for  X^,  X^,  and  y,  into  (28),  and  (40),  the 
optimal  lift  coefficient  CT  as  a  function  of  k  is  given  parametrically  by 


1 

.1  I  L  2 

CL=±2[8tan  (<r 2)/ct2]  /f(cr) 


(41a) 


K  - 


-1  ± 


,(l^)tan~  (a2  ) 


1 


i  i  I 
cr2 tan"  (a2) 


(41b) 


where  the  common  denominator  is 

8  tan  1(<Ja  ) 


f  (<*)  =  3± 


1 

a2 


a2  tan-1  a2  ) 


1_ 

a 


(41c) 


The  maximum  lift  coefficient  (+  sign  in  (41)  )  as  a  function  of  k  is  shown  in 
Fig.  2.  The  following  limiting  cases  are  of  special  interest. 


(ia)  K  0+  (a-*-  +  oo ) 

By  reducing  k  to  zero,  the  plate  is  constrained  to  be  nearly 
stretched  straight,  implying  that  T  and  p  =  O(r)  are  both  small.  Setting 
=  «  »  1,  we  deduce  from  (39),  and  (41b,  c),  that 

^  "  r  -  7 +°(€3  >  *  y  [i  +  ~  +  o(<2>]  (42) 

and 


€  =  (3«)*[l±  y  ( 3ti2  k  )*  +O0C*)]  . 

Finally,  the  optimal  lift  coefficient,  by  (41a),  is 

CL  =±  y  (SJKf'l  +J  (3TTZKf  +  0(jt*)]  . 

This  optimum  solution  shows  that  as  soon  as  K  increases  from  zero,  the 
plate  starts  to  bend  more  near  the  two  ends  and  carries  most  of  the  lift 
there,  since,  by  (36)  and  (42),  both  F  and  p  have  square-root  singularities 
at  |  =±  1.  As  the  arc-length  S  further  increases  for  fixed  chord  i  ,  k 
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and  «  become  greater,  \i  smaller,  so  that  the  singularities  at  the  ends 
then  become  weaker.  The  rate  of  increase  of  the  optimum  lift  with  respect 

-V4 

to  increasing  K,  dC^/d/c,  behaves  like  K  as  K  0  , 

(  ib)  0  <  a  «  1 . 

Another  interesting  limit  is  as  a  0+,  in  which  case  we  find 

M  =  TF  "  T  a  +  °{°*)]  »  Y  =±\/T[l  +  ^-  a+O^a2)]  (43a) 

K  =  \  (5  +  3^/2) [l  -  (4  ifi)  +0'a^]  .  (43b) 

Here,  the  upper  (or  lower)  sign  is  for  the  maximum  (or  minimum)  lift  co¬ 
efficient,  which  is  obtained  from  (41a,  c)  as 

CL=±  J  (S\fz  +  6)[l  -  ^  (5  +  3/2) +  0  (a2)]  .  (44) 

In  the  limit  as  a  —  0+,  the  maximum  lift  coefficient 

CT  =  0.918(1 -0.076a)  at  k  =  0. 1082(1  -0. 77ta;  (45) 

±j  max 

appears  to  be  very  respectable,  particularly  for  such  a  small  K .  On 
the  other  hand,  it  seems  rather  questionable  whether  the  minimum  lift  co¬ 
efficient 


CLmin='11-2°^1"°-924CT^  at  K  =  1 . 320(1 -1 . 624a )  (46) 


can  even  be  close  to  any  physically  realizable  situation  because  the  low 
pressure  underneath  the  plate  would  most  likely  cause  air  to  ventilate  the 
entire  lower  surface  of  the  plate. 


The  corresponding  solution  of  T  and  j3  can  be  immediately  deduced 
from  (36)  as 


ns) 


+ 


i-s 

T+f 


2  z 
"  T 


+  0(</) 


(47) 
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a 

F 


fc  (log ’] 


+  Ofa2) 


(48) 


Thus,  in  this  limit,  the  flow  angle  p(£)  has  a  logarithmic  singularity  at 
the  end  points  £  =  ±  1 ,  which  are  of  the  weakest  type  singularity  in  this 
linearized  theory.  This  particular  profile  of  the  plate  results  in  a  constant 
pressure  distribution  over  the  plate  except  for  a  higher  order  logarithmic 
singularity  at  the  two  edges.  The  over-all  features  of  the  solution  in  the 
present  limiting  case  therefore  indicate  that  this  is  the  most  favorable 
optimum  state  arrived  at  when  the  variation  of  parameter  /c  is  considered 
in  this  final  stage. 


Case  (ii)  -1  <a  <  0.  The  ranges  of  X  and  X  now  become  either 
-  1  2 

(X^O,  Xz>  0)  or  (Xt<  0,  X2<  0)  ,  (49) 

corresponding,  respectively,  to  the  maximum  and  minimum  lift.  In  terms 
of  the  parameter  v  defined  by 

i 

tanh  vir  =  { -  <j)2  (v  >  0)  ,  (50) 


the  solution  of  (11)  and  (32)  for  this  case  is  found  as 


r(£)  =  y  cosh  vit  cos 


i  1-6 

viog  yyf- 


P(£)  =  y  cosh  vir  coth  vtt  sin  |v  log  » 


(51) 

(52) 


where  the  coefficient  y  is  given  by  (35).  This  solution  can  be  immediately 
verified  by  making  use  of  the  formulae 


fl^sin  jvlog  ^ 

Hg  |cos  Jv  log  y^~ 


=  csch  vtt  -  coth  vir  cos 


vlog!+f  ’ 


=  coth  vtt  sin 


v  log 


!-e 

i+r 


(53) 

e i  <  i )  .  (54) 


The  above  inversion  formulae  can  be  derived  by  a  contour  integration  of 
(t-l)“'(t+l)  lv(t-£)  encircling  the  real  axis  from  t  =  -1  to  1  in  the 
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complex  t-plane. 

Although  both  of  the  above  T  and  p  remain  bounded  in 
-1  4  t  ^  1,  they  nevertheless  oscillate  infinitely  fast  as  the  end  points 
£  =±  1  ?re  approached,  and  hence  do  not  satisfy  the  inequality  condition 
(7),  which  is  imposed  on  physical  grounds.  For  this  reason  the  above 
solution  is  regarded  as  void  of  any  physical  significance,  and  hence  will 
not  be  further  pursued  here,  although  it  also  tends  to  the  limiting  solution 
(47),  (48)  as  ■-*•0.  However,  mentioning  of  this  case  may  serve  a  use¬ 
ful  purpose  to  point  out  that  when  a  numerical  iteration  method  is  employ¬ 
ed,  particularly  for  small  values  of  <j  ,  the  iterated  solutions  may  oscil¬ 
late  between  the  two  cases  (ib)  and  (ii)  and  the  success  of  such  procedure 
may  be  hindered  by  the  lack  of  convergence. 


5.  Discretized  Fourier  Series  Expansions 

We  next  consider  a  method  for  obtaining  approximate  solutions  to 
the  optimum  shape  problem  by  expanding  F{§)  and  P(£)  in  finite  Fourier 
series  with  the  coefficients  chosen  so  that  the  lift  is  maximized  under  the 
previously  mentioned  isoperimetric  constraints  of  fixed  chord  and  fixed 
wetted  arc-length.  This  Rayleigh-Ritz  method  of  solution  has  been  dis¬ 
cussed  previously  in  Part  I. 

Let  the  expansion  for  r(§ )  be  given  by 

N 

w-l  ansin(2n-l)6  ,  (55) 

n=l 

where  §  =  cos  0(0  ^  0  ^  ir)  and  the  an's  are  arbitrary  real  constants. 
This  choice  for  T  automatically  satisfies  the  end  conditions  (8).  From 
(11),  p(|)  is  given  by 

N 

P(§)  =  “^  ancos(2n-l)0  ,  (56) 

n=l 


where  we  have  changed  variables  (£  =  cos  0  ,  t  =  cos  <p)  and  made  use  of 
the  identity 


L  sin  rruy  sin <p 
J  cos  ip  -  cos 


sul<?T  d <p  =  -ir  cos m0 


o 
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Note  that  even  order  Fourier  components  have  been  omitted  from  the 
above  expansions  to  satisfy  the  symmetry  properties  (21 ). 

The  problem  of  maximizing  the  lift  for  fixed  arc-length  and  chord 

is  equivalent  to  minimizing  the  function  I  in  (22),  which,  by  (55)  and  (56), 

may  now  be  considered  to  be  an  ordinary  function  of  the  coefficients 

{a  }  and  the  factor  A.  To  minimize  I  we  set  the  partial  derivatives  of 
n 

I  with  respect  to  A  and  each  of  the  a^'s  to  zero,  giving  (since  A  does 
not  vanish) 


j  jr  it 

X  ~  e^cos  p  sin0d0  +  X  ^  e^sin0d0 


a  =  0 


^  -gO-  =  X  ^  e^sin[  p  +(2n-l  )0]  sin  0  d0 


r% 

+  X^  e^sin(2n-l  )0  sin0  d0  -  ^  6^  =0  n  =  1,2,  . 


•  ■  »  N  . 


(58(n)  ) 


The  (N+l ) -equations  above,  in  which  T  and  |3  are  given  by  (55)  and  (56), 
together  with  (27)  in  which  K  is  given,  determine  the  N  Fourier  co¬ 
efficients  (an)  and  the  two  Lagrange  multipliers  X^  ,  X^  .  Finally,  the 
lift  coefficient  is  given  by  (28).  Since  (57)  and  (58)  are  transcendental  in 
the  unknown  Fourier  coefficients  the  solutions  must  be  found  numerically. 

In  order  to  illustrate  this  method  we  consider  the  case  N  =  2. 
Equations  (58)  read  as  follows 


c  c 

X^J  e^sin(p +0)sin 0  d0  +  X^  J  e^sin^d©  =  j- 


(58(1)  ) 


X  \  e  sin(p +30)sin0  d0  +  X 


sin  30  sin  0  d0  =  0 


(58(2)  ) 


These  equations  are  to  be  solved  together  with  (27)  and  (57)  for  a  ,a  ,  X  , 

12  1 

and  X^.  The  Lagrange  multipliers  maybe  eliminated  from  (57),  (58(1)  ), 
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and  (58(2)  )  to  give  a  relation  between  a  and  a;  namely, 

1  2 

rX  X 

f(a^,  a^)  =  a^J  j  e^sin(p  +0)sin0  d0  J  e^sin  30  sin0  d0 


T 

-f 


e  sin({3 +30)sin  0  d6 
o  o 


.*  r  2  1 

j  e  8in20  d0 


-  ^  e^sin30  sin0  d0  ^  eicosp  sin0d0 


.7  w 

+  ^  e^8in(P  +30)8in  0  d.9  e^sin0  d0  =  0 


(59) 


A  second  relation  between  a  and  a  is  given  by  (27b)  in  whicl  k  is 
given  a  specific  value;  however,  it  is  somewhat  easier  to  do  the  inverse 
problem  in  which  k  is  determined  once  a  is  known.  Thus,  if  a^  is 
given,  a^  is  determined  by  (59),  and  K  is  then  fixed  by  (27b).  Finally, 
Cj^  may  be  found  from  (28),  (20),  and  (16)  (with  a  =0). 

The  curve  of  a  versus  l  satisfying  (59)  is  shown  in  Fig.  3  and 
the  lift  coefficient  for  the  two -term  Fourier  expansion  is  plotted  in  Fig.  2. 
As  a  0,  it  can  be  shown  by  expansion  of  (27b),  (28),  and  (59)  for 
I  aj  ,  faj  «  1,  that 

a  =  yy  a  +  0(a2  )  , 

2  17  i  i 

and 

CLmax  =  J\/t TT  ^  +  0('cVZ)  aS  K  ~*°  ■ 

Thus,  the  maximum  lift  coefficient  increases  more  slowly  (with  increas¬ 
ing  k )  for  the  two-term  Fourier  expansion  than  for  the  linearized  theory 
of  the  previous  section  (see  (42a)  ).  This  is  thought  to  be  due  to  the  pre¬ 
viously  mentioned  endpoint  singularities  which  are  present  in  the  linear¬ 
ized  theory. 

Actual  plate  shapes  for  the  case  N  =  2  are  shown  in  Fig.  4  for 
various  values  of  k  (note  change  of  vertical  scale).  These  are  found  by 
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numerically  evaluating  (15)  for  real  £.  The  factor  A  in  (15)  drops  out 
after  normalization  of  the  chord  to  unity.  Note  that  the  maximum  height 
of  the  plate  occurs  at  two  symmetrically  located  points  (at  the  X's  in 
Fig.  4)  for  smaller  values  of  k  (<,  0.030). 

The  cases  N  =  3,4,  .  .  .  ,  etc.,  could,  in  theory,  be  solved  as 
outlined  above  and  should  result  in  higher  and  higher  lift  coefficients  for 
a  given  ratio  of  arc-length  to  chord;  however,  the  calculational  difficulties 
involved  in  the  solution  of  the  system  of  equations  (57)  and  (58)  would 
surely  increase.  Fortunately,  the  numerical  experiments  in  Part  I 
indicate  that  just  a  few  terms  in  the  Fourier  expansions  are  needed  to 
give  results  which  are  within  one  per  cent  of  exact  solutions  in  those 
specific  examples. 
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Fig.  1  The  physic 


Optimum  lift  coefficient  C.  versus  k  =  S/l  -  1  for  t 
linearized  theory  and  the  tw6-term,  nonlinear  Fourier 

series  expansion. 


Optimum  planing  surface  shapes  for  the  two-term  Fourier 
expansion  (x's  denote  maximum  heights). 


THEORY  OF  OPTIMUM  SHAPES  IN  FREE -SURFACE  FLOWS 


PART  HI 

Minimum  Drag  Profiles  in  Infinite  Cavity  Flow 


by 

Arthur  K.  Whitney 

California  Ii  itute  of  Technology 
Pasadena,  California 

The  problem  considered  here  js  that  of  determining  the  shape  of  a  sym¬ 
metric  two-dimensional  plate  so  that  the  drag  of  this  plate  in  infinite 
cavity  flow  is  a  minimum.  The  flow  is  assumed  to  be  steady  and  irrota- 
tional;  effects  due  to  gravity  are  ignored.  With  the  aid  of  a  recently 
developed  theory  of  variational  calculus  involving  singular  integral  equa¬ 
tions,  the  drag  of  the  plate  is  minimized  under  the  constraints  that  the 
width  and  wetted  arc  length  of  the  plate  are  fixed.  The  extremization 
process  yields,  in  anology  with  the  classical  Euler  differential  equation, 
a  pair  of  coupled  nonlinear  singular  integral  equations.  Although  analytical 
and  numerical  attempts  to  solve  these  equations  prove  to  be  unsuccessful, 
analysis  of  the  equations  shows  that  the  optimal  plate  shapes  must  have 
blunt  noses.  Finally,  optimal  shapes  are  obtained  by  Fourier  series 
expansions  for  various  ratios  of  arc  length  to  plate  width. 


1 .  Introduction 

We  consider  the  two-dimensional  cavity  flow  of  an  incompressible 
fluid  past  a  plate  of  arbitrary  shape.  The  flow  far  upstream  is  uniform 
with  velocity  U,  pressure  p^,  and  density  p.  The  pressure  pc  inside 
the  cavity  is  assumed  to  be  a  constant,  so  that  by  Bernoulli's  law  the  fluid 
velocity  at  the  cavity  wall  is  a  constant  V,  where 


i  -Tirw«'>rr*'WFi  V’ 
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7  »V2  +  Pc  =  j  p'^ PB  • 

The  cavity  flow  may  be  characterized  by  the  non-dimensional  cavitation 
number 

“  =  <p«-  Pc*'  r  p,j!  • 

As  the  cavitation  number  decreases,  the  length  and  width  of  the  cavity 
grow  indefinitely  and  the  flow  approaches  the  "Helmholtz  flow",  in  which 
the  cavity  is  infinitely  long  and  the  cavity  pressure  equals  the  free  stream 
pressure  (ct  =  0)  so  that  the  cavity  is  maintained  far  downstream. 

The  specific  problem  considered  here  is  to  find  the  shape  of  a 
symmetric  plate  (see  Fig.  1)  of  given  wetted  arc  length  Sq  and  given 
width  y  ,  so  that  the  drag  of  this  plate  in  infinite  cavity  flow  (cr  =  0) 
is  a  minimum.  A  precise  definition  of  the  class  of  plates  under  consider¬ 
ation  will  be  given  in  the  next  section.  The  solution  of  this  problem  has 
obvious  applications  in  the  design  of  struts  or  other  two-dimensional  non¬ 
lifting  surface  which  may  operate  in  the  super -cavitating  range.  For 
large  Reynold's  number  flows,  the  viscous  effects  may  be  ignored  as  a 
first  approximation;  however,  corrections  due  to  viscous  drag  can  be 
calculated  once  the  potential  flow  is  kno  Finally,  although  optimal 
shapes  are  sought  for  the  case  ct  =  0,  these  shapes  should  be  approximate¬ 
ly  the  same  as  those  which  would  be  found  for  a  >  0.  This  can  be  seen  by 
considering  the  approximate  rule,  Cq(ct)  =  C^(0)(l+a),  which  relates  the 
drag  coefficient  at  CT  =  0  to  the  drag  coefficient  for  ct  >  0.  (Here, 

CT  <  I,  see,  e.g.  Gilbarg  I960).  Thus,  to  minimize  at  a  given 

ct  >  0,  we  could  just  as  well  minimize  0^(0).  It  should  be  noted  that 
this  rule  appears  to  be  a  good  approximation  only  for  blunt  bodies,  so  the 
above  argument  may  be  limited  to  the  case  when  the  arc  length  is  not 
much  larger  than  the  width  of  the  plate. 

Lavrentieff  (1938)  gave  the  solution  to  a  related  minimum 
drag  problem  (see,  e.g.  Gilbarg  (I960)  ),  that  of  finding  the  shape  of  a 
symmetric  plate  of  minimum  drag,  although,  in  Lavrentieff 's  problem, 
the  plate  is  confined  to  lie  within  a  rectangle  which  circumscribes  the 
nose  and  ends  of  the  plate.  If  the  nose  of  the  plate  is  at  (0,  0)  and  the 
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ends  of  the  plate  are  at  (xq,  ±  yQ/2),  the  solution  for  the  optimal  profile 
(see  Fig.  2)  is  found  to  consist  of  a  straight  section,  from  (0,  -h/2)  to 
(0,h/2),  and  the  free  streamlines  which  leave  the  ends  of  this  section  and 
go  on  to  pass  dirough  (xq,  ±  yQ/2).  The  length  of  the  flat  nose  section, 
h,  is  uniquely  determined  by  the  given  xQ/yo-  The  pressure  difference 
across  the  flat  portion  of  the  plate  is  the  only  contribution  to  the  drag 
since  the  fluid  pressure  equals  the  cavity  pressure  (p  =  pc)  on  the  free 
streamlined  sections  of  the  plate.  This  solution  was  obtained  by  the  use 
of  several  comparison  and  monotonicity  theorems  which  follow  from  the 
maximum  principle  for  harmonic  functions. 

The  present  work  was  originally  conceived  as  a  confirmation  of 
Lavrentieff's  solution,  but  it  was  hoped  that  this  could  be  done  by  using 
the  variational  calculus  method  introduced  earlier  in  Parts  I  and  II;  how¬ 
ever,  no  satisfactory  method  was  found  for  imposing  the  condition  that 
the  plate  be  confined  to  lie  within  the  rectangle.  On  the  other  hand,  if 
this  constraint  is  dropped,  one  can  easily  construct  a  sequence  of  plate 
shapes  which,  in  the  limit,  have  zero  drag,  neglecting  viscous  effects, 
of  course.  Such  a  sequence  is  illustrated  in  Fig.  3.  A  typical  plate 
consists  of  an  inverted  cap  of  width  h  and  length  plus  the  free  stream¬ 
lines  which  issue  from  the  ends  of  the  cap  and  go  on  to  pass  through  the 
corners  of  the  rectangle.  All  other  plate  shapes  in  this  sequence  are 
found  by  decreasing  h^  and  increasing  h^  in  such  a  way  that  the  free 
streamlines  still  pass  through  (xq,  ±  yQ/2).  As  h^ -*■  0,  the  flow  inside 

the  cap  becomes  a  dead  water  region  with  stagnation  pressure 

1  2  1  2 
Ps=  2  P^  +  P*.  so  that  the  drag  of  the  plate  is  just  ^  pU  h  ,  whicn  can 

be  made  very  small  by  choosing  h  sufficiently  small.  Note  that  the 

pressure  difference  across  the  back  face  of  the  cap  is  the  only  contrib’v 

tion  to  the  drag  of  the  plate. 

This  observation  led  the  author  to  consider  the  problem  described 
earlier.  It  was  thought  that  by  fixing  the  arc  length  of  the  plate,  shapes 
such  as  those  described  above  would  be  eliminated.  The  constraint  on 
the  length  of  the  plate  was  dropped  for  simplification. 
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2.  The  Problem  of  the  Symmetric  Cavitating  Plate 

The  class  of  flows  under  consideration  is  limited  to  infinite  cavity- 
flows  past  plates  P  with  the  following  properties  (see  Fig.  1): 

(a)  The  width  of  P  is  y  . 

(b)  The  arc  length  of  P  is  sq. 

(c)  P  has  a  continuous  slope  except  at  the  nose  wh  .re  the  vertex 
angle  is  2 a,  with  0 it. 

(d)  Let  S  and  S'  be  points  on  the  intervals  OA  and  OA', 
respectively.  The  pressure  p  on  S'OS  satisfies  p^Pc> 
while  on  SA  and  S'A',  p  =  pc;  that  is,  SA  and  S'A'  are 
free  streamlines. 

(e)  P  is  coincident  with  the  free  streamlines  which  issue  from 
S  and  S'. 

The  condition  p  ^  Pc  in  (d)  is  an  obvious  statement  of  the  fact 
that  the  vapor  pressure  (assximed  equal  pc)  is  the  minimum  pressure 
in  the  flow.  Sections  SA  and  S'.x'  are  included  since  free  streamlines 
have  already  been  shown  to  mak”  up  part  of  Lavrentieff 's  profiles  and 
similar  results  are  expected  for  the  present  problem.  It  is  more  con¬ 
venient  to  account  for  this  expectation  from  the  beginning  than  not.  Note, 
also,  that  (d)  imposes  no  midue  restrictions  on  the  problem,  since  the 
actual  locations  of  S  and  S'  are  not  known,  a  priori,  but  must  be  found 
as  part  of  the  optimization  process. 

Condition  (e)  guarantees  that  the  shape  ex  an  optimal  profile  will 
be  unique,  since  the  plate  surfaces  SA  and  S'A1  could  obviously  be 
moved  inward  toward  the  x-axis  without  changing  the  overall  forces  on 
the  plate.  This  condition  simply  provides  a  one-to-one  relationship  be¬ 
tween  plate  shapes  and  the  resultant  flow  patterns.  Note  that  there  are 
many  ways  of  doing  this;  the  plate  surface  could,  for  example,  be  required 
to  run  along  the  straight  lines  joining  S,  A  and  S',  A1. 

By  proper  choice  of  origin  and  magnification,  the  complex  potential 
plane  f  =  cp  +  i»J>  is  mapped  to  the  upper  half  t,  =  %  +  ir|  plane  (see  Fig, 

4)  by 
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f  =  j  All?,2  ,  (1 ) 

where  q>  is  the  velocity  potential,  is  the  stream  function,  and  A  is  a 
real,  positive  constant  which  is  chosen  so  that  S'OS  maps  to  |  £  |  ■$  1 . 

The  sections  SA  and  S'A*  map  to  1  <  |£|<  c,  where  c^l  (equality 
only  if  S  =  A,  S'  =  A1 );  the  remaining  sections  of  the  free  streamlines 
lie  on  c  ^  |  %  |  <  oo . 

We  next  introduce  the  complex  velocity 

w  =  df/dz  =  u  -  iv  =  qe  ^  (2) 

and  the  logarithmic  hodograph  variable 

«(?,)  =  log(U/w)  =  log(U/q)  +  id 

=  T(|,q)  +  i0(|,q)  ,  (3) 

where  u  and  v  are  the  x  and  y  components  of  the  fluid  velocity, 
q  =  (u2+  v2)2  is  the  speed,  and  9  is  the  direction  of  the  flow  with  respect 
to  the  positive  x-axis.  With  these  definitions,  Bernoulli's  law  may  be 
written  as 

P  -  Pc  =  \  P(U2  -q2)  -  \  pU2  (1  -  e'2  T)  ,  (4) 

where  p  is  the  pressure  at  any  point  in  the  flow.  On  the  free  streamlines 
SI  and  S'l1,  p  =  p£;  therefore,  by  (4), 

t(£,0+)  =  0  ,  |el>l  •  (5) 

If  the  boundary  value  of  w  on  S'OS  is  denoted  by 

W(e+iO)  =  r(e,o+)+ie(e,o+)  =  r(e)+ip(e)  ,  UUi  .  (6) 

where  F  and  p  are  real,  then,  since  pc  is  the  minimum  pressure  in 
the  flow,  the  inequality 


r(£)^o 


UU  i 


(?) 
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follows  by  (4),  and,  because  the  pressure  is  continuous  at  S  and  S', 

r  (±1)  =  0  .  (8) 

The  hodograph  variable  may  be  split  into  two  parts, 

«(£)  =  ^^(4)  +  «(£),  where  accounts  for  the  singular  behavior  of  w 
at  the  stagnation  point,  £,  =  0,  and  u(=r  +i(3  ,  for  |£|  <  1) is  analytic  in 
the  entire  upper  half  plane  including  the  real  axis.  It  can  be  shown,  in 
fact,  that 

«(l)  =  ^  log  {[  (C2  -i)^  +  i]/U  +  «(£>  ,  (9) 

where  (t,2  -l)2  is  taken  to  be  cut  along  the  real  axis  |£|  ^  1  and  is 
positive  for  L,  =  £  >  1;  the  logarithm  function  is  defined  to  be  that  branch 
which  is  real  for  a  real,  positive  argument,  with  a  cut  along  the  nega¬ 
tive  real  axis  of  the  argument.  Letting  t,  -»£  +  iO,  |f||  ^  1,  in  (9)  and 
comparing  this  ",'ith  (6)  we  have 


r<e)  =  ~  iog{[i+(i-e2)2]/le|}  +rte) 

(10a) 

IMS)  =  a  sgn£  +  p^)  , 

(10b) 

so  that  the  proper  logarithmic  singularity  of  T  =  Rew  =  log(U/q)  and  the 
proper  jump  in  flow  angle,  (3(0+)  -  (3  (0-)  =  2a,  are  exhibited  at  t,  =  0. 
Since  the  real  part  of  the  first  term  of  (9)  vanishes  for  £,  real,  |£  |  >  1, 
we  have,  by  (5), 

r(|)  =  Re{w(£+i0)}  =  0  ,  lei  >1  •  (ID 

Furthermore,  the  flow  approaches  that  of  the  free  stream  (q  =  U,  0  =  0) 
at  large  distances  from  the  plate;  therefore,  w  0  as  |^|  -*■00,  by  (3), 
and 

«(&)-0  ,  as  IU  -*«>  ,  (12) 

by  (9).  Finally,  r  in  (10a)  must  be  chosen  so  that  inequality  (7)  is 
satisfied  and  (8)  and  (10a)  imply 

r  (±1)  =  o  . 


(13) 
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The  function  co  may  be  determined  for  a  given  plate  shape,  how¬ 
ever,  because  the  solution  involves  complicated  nonlinear  functional 
equations,  it  is  easier  to  deal  with  the  so-called  "inverse  problem",  in 
which  the  plate  shape  is  determined  by  giving  either  T  (£)  or  (M£). 

This  information,  together  with  (11)  and  (12),  determines  gj^  uniquely. 

We  now  suppose  that  T  =  Re  go  is  given  which  satisfies  (13)  and  is  Holder 

X  1  1 

continuous'  on  |  £  |  1 .  The  Dirichlet  problem  for  the  determination  of 


Re{c^  (£+i0)}  = 


F(g)  for  UUl 


!§l  >  1  , 


together  with  (12),  has  the  solution 

(  .'■I  r(t)dt 

TT-  •  (,5) 

-1 

which  may  be  verified  by  letting  £  -*•§  +  i0(  1 1  j  ^  1 )  in  (1 5)  and  using 

Plemelj's  formula  (e.g.  ,  see  Muskhelishvili  (1953),  §17).  The  imaginary 

part  of  go  on  S'OS  is  found  to  be 
l 

,  r1  r(t)dt 

P1(S)  =  -  *t_|  ,  IsUi  .  (16) 

-l 


where  C^.  denotes  the  Cauchy  principal  value.  Holder  continuity  of  (3 
follows  from  the  assumptions  which  have  been  made  on  T  .  (Note  that  if 
T  does  not  satisfy  (13),  but  approaches  a  nonzero  value  at  an  endpoint, 
then  p  ,  as  given  by  (16),  will  have  a  logarithmic  singularity  at  that 
endpoint.  )  If,  instead  of  T,  p  is  given,  the  solution  for  to  involves  a 
Riemann-Hilbert  problem;  however,  this  solution  will  be  omitted  since  it 
will  not  be  needed. 

By  (2)  and  (3),  the  physical  plane  is  obtained  by  integrating 

dz  =  i  eWdf  =  Ae^^t,  dt,  .  (17) 


TF(£)  is  said  to  be  Holder  continuous  on  [  -1,1]  if,  for  any  two  points 

€,.£,€[ -1’1!’  |r(£)-nu|  <B|g  -g  I**,  with  B>0  and  0<JX<1. 

1*  1112  12 
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Thus,  the  plate  shape  is  given  parametrically  by 

£ 

z(l)  =x(|)  +  iy(|)  =  e“H  dt,  ,  UUc 


and,  since  the  plate  is  symmetric,  the  width  is  given  by 


yo=ImAj 


e“^  dt, 


(18) 


-c 


It  is  convenient  for  subsequent  analysis  to  convert  this  expression 
for  the  width  to  an  integral  from  -  1  to  il.  To  do  this,  we  first  continue 
w(4)  into  the  lower  half  plane  by  to(^)  =  -to  KV,  so  that 
w(£±iO)  =±r(£)+ip(£)  for  |||<tl  and,  by  (5),  to(£+iO)  =  to(£  -iO)  for 
|$  |  >1.  Next,  we  consider  the  function  which  appears  in  the 

integrand  in  (18).  This  function  is  uniquely  determined  by  the  jump  in  its 
value  (due  to  the  discontinuity  of  to)  across  the  cut  |£|  ^  1  and  by  its 
expansion  for  large  1 |  (see  Muskhelishvili  (1953),  §78).  In  fact,  it  can 
be  shown  that 


=  - - t‘-T"  ^  j  r(t)dt  +  6.  ,  (19) 

-1  -1 

in  which  the  first  integral  exhibits  the  correct  discontinuity  across  the  cut 
j  £  |  1  and  the  last  two  terms  are  required  by  the  expansion  of  the  left  side  of 

(19)for  |^|»1  (see  Whitney  (1969) ).  In  this  expansion,  we  use  the  integral 
representation  for  w,  r'  1 

-1 


n  te^  ^sinh  F  (t)dt 


(20) 


which  is  found  by  exactly  the  same  procedure  as  that  used  in  determining 
w  (&)  (see  (14)  and  (15)  ).  By  substituting  (19)  into  (18),  and  noting  that 
p  (-£)  =  -  P(£),  we  obtain 

i 


y  =  — 
7  o  tr 


1 


J  t  sinhr(t)sin  p  (t)dt  +  ^  t  sinhr(t)cos  p  (t)log{(c-t)/(c+t)}dt 


-1 


-1 

1 

+  2c  ^  r (t)dt 

-1 


(21) 


64 


An  element  of  plate  arc  length  ds  is  found  from  (17),  (5),  and 
(6),  to  be 


'Aer(e)|£id£  for  |$Ul 


ds  =  dz  =< 


L|£|d£ 


for  iei  >  i  . 


Thus,  the  total  arc  length  of  the  plate  is  given  by 
c  1 


s 


=  y  ds  =  aJ*(cz -1 )  +  ^  er^|t|dt 
-c  *"  -1 


(22) 


The  complex  force  acting  on  an  element  of  the  plate  dz  is  given 


by 


dF  =  (p  -  pc)(-idz)  , 


or,  by  (4 ), (5),  (6),  and  (17), 


dF  =< 


j  PU2A(l-e"Zr)er+lp£  d| 

o  ,  le I  >  i 


ill  ^  i 


By  integrating  this  expression,  we  obtain  the  drag  D  and  lift  L  acting 
on  the  plate  as 

1 

D  +  iL  pUzA  J  t(l-e'2r)er+lpdt  , 

-I 

so  that,  again  noting  the  asymmetry  of  (3  (£ ), 

1 

D  =  pUzA^  t  sinhr(t)sin|3  (t)dt  ,  (23) 

-1 

and  the  lift  vanishes,  as  we  should  expect  for  a  symmetric  shape  at  zero 
angle  of  incidence.  An  alternate  expression  for  the  integral  which  ap¬ 
pears  in  (23),  and  also  in  the  expression  for  the  width  (21),  may  be  ob¬ 
tained  by  substituting  (20)  in  the  left  side  of  (19)  and  expanding  this  equa¬ 
tion  for  large  |  £  j .  By  matching  coefficients  which  multiply  like  power 
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of  4  in  this  expansion,  we  obtain  the  identity 


1 

^  t  sinh  r  (t)sin  0  (t)dt 
-1 

- 1 

from  the  terms  multiplying  4 


2 

r(t)dt) 


(24) 


Since  p  and  U  are  kept  constant  in  the  minimization  of  the  drag, 
it  is  convenient  to  give  the  drag  the  dimensions  of  length  by  setting 
D  =  D /  £-  pU2.  By  (23)  and  (24)  we  have 

2 

r(t)dt)  .  (25) 

-i 


3.  Statement  of  the  Minimum  Drag  Problem 

In  the  previous  section  we  have  shown  that  the  problem  of  minimiz¬ 
ing  the  drag  of  a  symmetric  profile  of  given  width  and  arc  length  reduces 

to  finding  functions  T  (4)  and  (3  (4)  and  constants  A,  c,  and  a,  so 
#  li 

that  D  in  (25)  is  a  minimum  subject  to  the  constraints  (21)  and  (22),  in 
which  yQ  and  sq  are  given  fixed  quantities.  The  functions  r(|)  and 
0(4)  in  (21),  (22),  and  (25),  are  related  to  r(4),  0^),  and  the  half 
vertex  angle  of  the  plate  a  by  equations  (10);  furthermore,  0(1)  is 
related  to  r(4)  by  (16)  in  which  T  satisfies  the  end  conditions  (13)  and 
is  chosen  so  that  T  in  (10a)  is  positive.  Equivalently,  we  may  state 
this  problem  in  terms  of  r(4),  0  (4),  A,  and  c,  and  omit  further  refer¬ 
ence  to  r(4),  0^(4),  and  or,  however,  the  discontinuous  behavior  of  T 
and  0  ,  as  exhibited  in  (10),  should  still  be  remembered.  By  letting 
4  -*4  +  iO  (|4|  ^  1)  in  (20),  we  obtain  the  identity 

\ 

0(4)  =  -  I<£  =  -He[r(t)]  (26) 

-1 

which  may  also  be  verified  by  (10)  and  (16).  In  the  above,  the  symbol 
H|  denotes  the  finite  Hilbert  transform.  Finally,  in  (21),  (22),  and  (25), 
the  factor  A  is  a  real,  positive  constant  and  the  parameter  c,  which 
determines  the  location  of  the  endpoints  of  the  plate  in  the  4 -plane, 
satisfies  1  c  <  oo. 
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This  problem  is  equivalent  to  that  of  finding  a  pair  of  extremal 
arcs,  r{£)  and  p(£),  which  satisfy  (26)  and  minimize  the  functioned 


l[r(£),p(£);  A,  c]  =  D  -  XiSo  -  X2yo 
1  1 
=  a]  f(r(t),p(t),t;  xifx2,c)dt 
-1 

/.I  l2 


+  ^  (!  -Xz/2tt)|J  r(t)dt 


where,  by  (21),  (22),  (24)  and  (25), 


f(ni).p (£),£;  x t> x 2> c )  =  -x£  (c2 -1)  +  er(e)|el] 

X 

-  [I  suih T (£ )cos  p  (£ )log{c -£ )/{c+£ )}  +  2cT(£)] 


(27a) 


(27b) 


In  the  above,  the  integral  identity  (24)  has  been  used  in  the  expression  (21) 
for  yQ  and  X^,  X2  are  Lagrange  multipliers. 

The  geneial  variational  problem  of  this  type  has  been  investigated 
earlier  in  Parts  I  and  II.  For  the  present  problem  the  method  of  solution 
will  follow  the  same  approach  with  a  few  modifications.  Let  the  set 
{r(£ ),  p  (£),  A,  c}  denote  the  optimal  solution  and  let  {T(£ ),  p(£ ),  A,"c}  be 
an  arbitrary  neighboring  set  which  also  satisfy  (21),  (22),  and  (26).  The 
differences  6T(£)  =  T  -  T,  6p  (|)  =  p  -  p  ,  6A  =  A-  A,  6c  =  c  -  c,  form 
a  set  of  small  variations,  where  by  (26), 

6p<!)  =  -Hgtfint)]  .  (28) 

The  variation  of  the  functional  I,  about  the  optimal  value,  due 
to  the  variations  {6T,  6p  ,  6A,  6c}  is  given  by 


Al=l[r+6r,  p  +6p  ;  A+6 A,  c+6c]  -l[r,p;A,c] 


Expansion  of  the  above  expression  for  small  |  6 F  |  ,  |  6 p  |  ,  |  6A|  ,  and 
|  6c  j  ,  yields 
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Al  -  61  +  62I  +  63I  +  . 


where,  by  (27),  the  first  variation  61  and  second  variation  62I  are 
given  by 

6 1  =  eAjjftr.p.t;  X^X^cJdt  +  I  (l  -X2/2tt)J2| 


+  A 6c  »  f  (r,  p  ,  t;  X  ,  X  ,  c)dt 

v  ^  12 

+  A  jjff^r.P.t;  XifX,,c)+  |  (l-X2/2ir)j]6r(t) 


+  fp  (F,  p  ,  t;  XifX2,c)6p(t)>dt 


62I  =  26A^|[fr+  |  (l-X2/2ir)j]6r+  fp6p|dt 

+  26  A6c  y  fcdt  +  A(6c)2  j  f.cdt 
(* 

+  Aj  {frr(6r)2+2fr(36r6p  +fpp(6p)2>dt 
+  ^  (l-X2/2ir)|j  5,rdt 


In  the  above,  all  integrals  are  from  -1  to  +1,  subindices  denote  partial 
differentiation,  and  J  is  given  by 


r% 

j  =  j[r<t)]=^  r(t)dt 


For  I  to  be  a  minimum,  we  must  have  61  =  0  and  62I  >  0  for 
arbitrary  6r(£),  6p(£),  6A,  and  6c.  61  vanishes  if  the  terms  multiply¬ 
ing  6 A  and  6c  vanish  separately.  The  first  line  in  (29)  is  zero  if 


f(F(t),p(t),t;  X j, X 2> c )dt  +  i-  (1-X2/2tt)J2  =  0  , 


(32a) 
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or,  by  (27a)  and  (31), 


D*  =  X.  s  +  X.  V  ,  (32b) 

l  o  t  o 

sjc 

in  which  the  functional  forms  of  D  ,  s  ,  and  y  are  given  by  (25),  (22), 

o  o 

and  (21),  respectively.  The  second  line  of  (29)  is  zero  if,  by  (27b), 


y  fc(F(t),  (3  (t),  t;  X^X^cJdt 


=  2X 


_inhr(_t)cos_p_(t)_  +  r(t)l  dt  =  0 


-t2 


(33) 


Finally,  for  the  last  integral  of  (29),  we  substitute  (28)  for  6(3  and  then 
change  the  order  of  integration,  giving 

|  (1  -W2Tr)j]  6T (t)  +  fp  6(3  (t|  dt 

fr+  |  (1  -\  j2ir)J  +  H^[fp]|6r(t)dt  =  0  .  (34) 

; 

Now,  since  6T(£)  is  arbitrary,  we  obtain  the  nonlinear  singular  integral 
equation 


f^r.p.l;  X^.c)  +  H^[fp(r,p,t;  X^.c)] 


r(t)dt 


(35a) 


where,  by  (27b) 


X 

— •  [|  coshr(£)cos  (3  (!)log{(c-£)/(c+!)}  +  2cl 


X 

fp  =  |  sinhr(|)sin(3(£)log{(c-$)/(c+S)} 


(35b) 


This  integral  equation,  which  contains  X  ,X  ,  and  c,  as  parameters, 

1  2 

is  to  be  solved  together  with  the  linear  integral  equation  (26)  for  the 
extremal  arcs  T(|)  and  (3(|).  The  Lagrange  multipliers  X^X^  are 


I 
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determined  by  equations  (32)  and  (33).  The  final  parameter  c  is  most 
conveniently  determined  by  the  nondimens ional  ratio  sQfvo>  which,  of 
course,  is  a  known  constant.  Finally,  the  optimal  drag  coefficient  (based 
on  plate  width)  is  given  by 


CD  *  D/  j  pU*yo  =  D*/y0  . 


(36) 


This  optimal  drag  coefficient  is  a  minimum  if  the  second  variation 
62I  >  0,  which  reduces  to 

62I  =  Atfcf^f^dt  +  A  J{frr(6rf  +2frp6^6p  +  fp  p  (6p  )r  }dt 

+  (1 -X^/2ir)|y  6T  dtj  >  0  , 


(37) 


by  (30),  (32a),  (33),  and  (34).  Now,  the  last  term  in  (37)  is  due  to  the 
second  variation  of 

2 


-i 

which,  by  (24),  equals  the  second  variation  of 

1 

2A(1  -\J2tt)  ^  t  sinhr(t)sin  p  (t)dt  , 

-1 

so  that  (37)  may  also  be  written  as 

A(6c)Z  Jfccdt  +  Ai  (grrl6r)Z+zZrp6r6P  +gpp(6P)Z}dt>  o  ,(38) 


where  g  =  f  +  2A(1  -W2u)£  sinhr(£  )sinp  (£ ).  Since  A  >  0,  this  inequality 
holds  if  both  integrals  are  positive.  The  first  term  is  positive  when,  by 
(27), 


f  f  dt  =  -2X  +  \  t2sinhr(t)cosp(t)dt  >  Q 

J  CC  1  *  (c2-t2)2 


(39) 


In  Part  I,  it  was  shown  that  a  necessary  condition  for  the  second  integral 
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in  (38)  to  De  positive  is  that 

grr  +  gp p  =  I  >  0  *  for  Ul<!  » 

or  simply, 

*■  <  0  .  (40) 

since  F  is  real.  Once  an  optimal  solution  has  been  found,  conditions  (39) 
and  (40)  must  be  checked  to  determine  if  the  solution  is  an  actual  minimum. 

The  singular  integral  equation  (35)  may  be  reduced  to  an  integral 
equation  with  a  regular  kernel  by  using  the  identity 

£  coshr(£)cos  p  (£)  =  H^[  tsinhr(t)sinp  (t)]  +  £  ,  (41) 

which  follows  ^y  averaging  the  equations  which,  in  turn,  are  a  result  of 
letting  £  -*-£±i0(|£|  <  I)  in  (19)  and  noting  the  limits  u(£±i0)  =  ±r(£  )+ip  (£). 
The  substitution  of  £  coshT  cos  (3  ,  as  given  by  (41),  into  (35)  yields 

X.  \  t  sinhr(t)sin(3  (t)K(t,£ ;c)dt  -  X  er^  |£  | 

Z*J  1 

=  -i-  (2c  +  6 log{ (c -I )/ (c+g )>  -  |  (1-^/2*)^  T(t)dt  ,  (42a) 

-1 

where 

K(t.£;c)  =  -L  log{(c-t)(c+£)/(c+t)(c-£)}/(t-£)  .  (42b) 

TC 

It  is  now  possible  to  show  that  if  an  optimal  plate  shape  exists,  it 
must  have  a  blunt  nose;  i.  e.  ,  a  =  ir/2.  To  show  this,  we  first  note  that, 

•p  \ 

by  (42),  e  ^|£l  possesses  a  regular  series  expansion  about  £  =  0  of 
the  form 

er(e)  III  =  c  +  c£2  +  .  .  .  ,  (43) 

U  1 

where  the  explicit  forms  of  c  ,  c  ,  ,  etc.  ,  are  easily  found  but  will 

be  omitted  here.  On  the  other  hand,  from  (10a),  in  which  T(£)  is  con¬ 
tinuous  , 
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er(6>  lei  =  [{i+(i-e¥}/lel]2“'/’ieUI‘te) 

=  piw,I5I'‘w’  e^’  +  odel21'-0,7’1) 

as  |£  |  —■  0.  This  expansion  agrees  with  (43)  only  if  a  =  ir/2. 

Since  there  are  no  known  analytical  methods  for  solving  the  system 
of  integral  equations,  (26),  (42),  our  only  recourse  is  in  solution  by  numerical 
methods.  These  attempts  have  also  failed;  however,  in  order  to  illustrate 
some  of  the  many  difficulties  which  beset  such  procedures,  we  briefly  men¬ 
tion  one  of  the  schemes  that  has  been  tried.  First,  the  integrals  in  (26), 

(32),  (33),  and  (42),  are  approximated  by  numerical  quadratures  which  in¬ 
volve  the  values  of  F(£)  and  (3  (£ )  at  N  points  {£.}  from  -1  to  +1 .  An 
initial  guess  of  (r(£.)}  is  made  and  the  set  {p(§.)}  is  then  given  by  (26). 
Next,  the  Lagrange  multipliers,  X  ,\  ,  are  found  from  (32)  and  (33),  in 
which  the  current  values  of  T  and  p  are  used.  Finally,  new  values  of 
{F^)}  are  calculated  by  solving  for  F(£)  in  (42)  and  the  process  is 
repeated,  hopefully,  until  the  iteration  converge.  These  calculations  are 
done  for  arbitrary  values  of  c  ^  1,  with  different  c's  corresponding  to 

different  ratios  s  / y  . 

o  o 

As  mentioned  above,  this  method  and  many  others  like  it  do  not 
work.  Among  the  more  disagreeable  features  that  are  encountered  are  the 
following:  (i)  the  iterations  do  not  always  converge;  (ii)  the  values  of 
r(§^),  as  given  by  the  solution  of  (42)  for  T(£),  are  not  always  positive, 
so  that  (7)  is  violated;  and,  (iii)  the  value  of  F  at  £  =  ±  1  is  not  zero; 
so  that  p  in  (26)  becomes  large  as  £  1  and  integration  of  terms  in¬ 

volving  sinp  and  cosp  by  numerical  quadrature  fails.  Corrective  steps 
were  taken,  such  as  enforcing  T(±l)  =  0  at  each  iteration  step,  and  in¬ 
creasing  N;  however,  these  measures  did  not  help. 

Note  that  P(|)  in  (26)  has  a  logarithmic  singularity  at  |  =  ±  1  if 
r(±l)  =  const  $  0,  and  p  has  a  higher  order  (branch  type)  singularity 
if  r  itself  is  singular  at  the  endpoints  (e.g.  ,  see  Muskhelishvili  (1953), 

§29).  By  studying  the  linear  case,  in  which  the  functionals  are  quadratic 
in  r  and  p  (see  Parts  I  and  II),  it  was  deduced  that  the  endpoint 
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condition  (7)  cannot,  in  general,  be  satisfied.  If  this  conclusion  is  also 
true  for  the  present  nonlinear  case,  it  is  likely  that  a  solution  to  (26)  and 
(42)  will  have  no  direct  physical  relevance  since  T  and  (3  will  be  singular 
at  £  =±  1;  nevertheless,  such  a  solution  (if  one  exists)  would  provide  an 
absolute  lower  bound  for  the  drag  which  could  then  be  used  in  judging  the 
"degree  of  optimality"  of  results  obtained  by  other  (approximate)  methods, 
such  as  that  to  be  presented  in  the  next  section. 


5 .  Solution  by  Expansions  in  Finite  Fourier  Series 

We  now  investigate  a  method  for  obtaining  approximate  optimal  solu¬ 
tions  by  the  expansion  of  r  (£ )  and  (3^(£)  in  Fourier  series  in  which  the 
constant  coefficients  are  chosen  so  that  the  drag  is  minimized,  subject  to 
the  isoperimetric  constraints  mentioned  previously.  Let  the  expansion  for 
r  be  given  by 


N 


^  =  "  2/  ansin'2n-1  • 


(44) 


n=l 


where  £  =  cos  d  (0<  6  ^  tt  .  This  6  is  not  to  be  confused  with  the  flow  angle). 
From  the  identity 

TT 

<£  sinmf  sinv%  =  -it  cos  m6  . 

J  cos  <p  -  cos  o 


we  see  that  (16)  is  satisfied,  term  by  term,  if 

N 


Pt(£)  =  anCos(2n-l)e  . 


n=l 


(45) 


Note  that  T  in  (44)  and  (3  in  (45)  have  the  correct  symmetry  properties 
and  that  (13)  is  automatically  satisfied.  By  setting  a  -  ir/2  in  (9)  (this  is 
dictated  by  the  results  of  the  previous  section)  we  obtain 

«(U  =  l°g{[(42-l)^+i]/U  +«(&)  . 

Therefore,  by  (17),  the  width  of  the  plate  is  given  by 
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f  “.(U  ,  i 

yo=-2AImj  el  [  (t?  -1  )2  +  i]  dt,  . 


This  integral  is  most  easily  evaluated  by  the  change  of  variables 

Z,  =  (v  +  v  )/2  ,  (47) 

which  maps  the  upper  half  t,  -plane  into  the  half  circle,  |  vj  1 ,  Im  v  ^  0 
(see  Figs.  4  and 5  ).  The  inverse  transform  v  =  -  t,  +(f,2  -l)2  is  chosen  so 
that  the  point  at  infinity  in  the  l,  -plane  maps  to  v  =  0;  the  endpoints  of  the 
plate,  Z,  =±  c,  maps  to  v  =  K,  where 

K  =  c  -  (c2  -  1)*  .  (48) 

It  is  readily  verified  that  u  as  a  function  of  v  is  given  by 


IN 

,i(v)  =  i£anv2n  ‘=ifl(v)  , 


*iv 

since  on  S'OS,  i>  =  ~e  ,  where  0^  0  <  it;  thus,  the  real  and  imaginary 
parts  of  (49)  agree  with  (44)  and  (45),  respectively.  From  (47),  (48),  and 
(49),  the  expression  for  the  width  becomes 


A 

yo  =  I 


Im  ^  e"  ^ v\ v+2i-2 v  -2iv  +v  )dv  , 


which  is  evaluated  by  taking  the  path  of  integration  L£  in  Fig.  5  In  the 
limit  e  -*■  0,  it  can  be  shown  (see  Whitney  (1969)  )  that 

r  * 

yo  =  j  y  \  sinQ(t)[2t  -  {2  -fi'(t>}2  It]  dt 


=  j  \  sinfi(t)[2t  -  {2-n'(t>}2  It)  dt 

i  ° 

+  j  cosfi(t)[  2  +  fi"(t)/2t]dt  +  J  (2-a  f 


+  [2//(  -  S2'(k)/2k]  cosB(k)  -  (l/2/cZ)sinf2(/c) 
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where  12' (t)  =  dft/dt,  fr"(t)  =  d2fi/dt?.  By  (10a),  (44),  and  (48),  the  expression 
for  the  arc  length  (22)  becomes 


.1  2 

(K  +  K  )  -1 


a 

n 


sin(2n-l  )0< 


(1+sin  0)sin  6  d0 


(51) 


Finally,  from  (10a),  (25),  and  (44),  the  drag  is  found  to  depend  only  on  the 
first  of  the  Fourier  coefficients, 

D*=  ^  (2-a)2  ,  (52) 

due  to  the  orthogonality  of  the  set  {sin(2n-l)0}  on  [  0,tt]  . 

The  optimization  problem  reduces  to  minimizing  D  in  (52),  sub¬ 
ject  to  the  constraints  (50)  and  (51),  over  the  (N+2)  -  dimensional  space 
(A.K.a^,  a^,  .  .  .  ,  a^).  For  general  values  of  N>1,  this  problem  must 
be  done  numerically;  however,  if  N  =  1  the  integrals  in  (50)  and  (51)  may 
be  evaluated  in  terms  of  special  functions.  Note  that  by  (10b)  and  (45), 
p(£)  =  y  sgn£  -  a^cosfl,  for  the  case  N  =  1,  so  the  plate  section  S'OS 
is  convex  or  concave  when  viewed  from  the  approaching  flow  as  a^  is 
positive  or  negative.  This  section  is  a  flat  plate,  corresponding  to  the 
Lavrentieff  profile  (see  Section  1 ),  if  a^  =  0. 

With  N  =  1  in  (49),  S2(t)  =  a^t,  fi'(t)  =  a^,  and  ft"(t)  =  0,  so  that 
(50)  becomes 

r  * 

yQ  =  j-  ^  sin(a^t){2t  -  (2-a^)2  /t}dt 

K  O 

+  2  J  cos(at)dt  +  £  (2  -  a  f 
o 

+  (2 /k  -  a/2jt)cos(aic)  -  (l/2/c*)sin(aK) 

A  r  -2  -l  -2 

=  ^*[{2a  +  4a^  -K  )sin(a  k) 

+  (4k"  1 -a«" 1 -2/ca  *)cos(a  k) 
i  l  i 

+  (2-a  f{v!l  -  &(a*)}  , 


(53) 
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where  Si  is  the  sine  integral, 

aw,]1  . 

o 

From  (51),  the  arc  length  is  given  by 

r.  2  p17  -a  sin© 

so  =  A  ( K+K~  )  -1  +  \  e  1  (l+sin©)sin<?d© 

o 

=  x  [kz+K~2  +  6  +  4tt(1  +a  _1  ){L«(a  )  -  I(a  )}  (54) 

4  11111 

-  4ir{Lo(a)  -  IQ(a)}]  , 

where  L.  and  L  are  the  modified  Struve  and  Bessel  functions,  respective¬ 
ly  ( e.g . ,  see  Abramowitz  and  Stegun  (1964)  ). 

The  problem  of  finding  the  optimal  plate  shape  from  the  class  of 
plates  with  N  =  1  in  (44),  (45)  is  equivalent  to  extremizing 

>jc 

I(A,  K,  a^)  =  D  (A,  a^)  -  X^A.K.a^ 


-\^o(A,K,a)  , 

& 

where,  as  before,  X^  and  X^  are  unknown  Lagrange  multipliers  and  D  , 

sq,  and  yQ  are  given  by  (52),  (54),  and  (53),  respectively.  If  I  is 

extremal;  the  three  partial  derivatives  IA,  I  ,  and  I  ,  must  vanish. 

A  K  a 

i 

This  gives  three  relations  among  the  quantities  A, fc,a  ,  X  ,  and  X^.  By 
eliminating  X  and  X^  from  these  three  equations  we  have 


AzA(a  ,/c)  =  D.[  s  y  -s  y  ] 
'  i  A  oa  7o/c  ok1  oa 

l  l 


■^a  ^  soAy ok  So/cyoA^  ^  ’ 

l 


where  the  partial  derivatives  may  be  found  from  (52),  (5*),  and  (54).  Since 
A  >  0,  we  must  have  A(a  ,/c)  =  0.  Let  the  solution  of  this  equation  be  given 
by  a  =  f(/c).  For  k  ~  l(c  ~1)  it  can  be  shown,  by  rather  long  and  tedious 
expansions  of  (52),  (53),  (54),  and  (55),  that 
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=  f(jc)  =  {8/(3ir+1 6)}(1  -*c)2  -  {24tr/(3ir+16)2  }(1  -K  f  +  0(1  -K  J4  .  (56) 

The  general  solution,  plotted  in  Fig.  6,  is  found  by  fixing  k  at  various 
values  between  0  and  1  and  numerically  solving  for  from  A(a  ,k)  ~  0 
in  (55).  A.s  K  0  (c  a  is  found  to  be  the  root  of  the  transcendental 

equation 

(2+a^)a  {Lo(a)  -  I^a^)}  -  (a2 +2a^ +4){Li(aj)  -  I^a^)} 

=  I  (3-2a)a2/(2-a)  . 
ir  ii  l 

This  root  is  given  by  a  C.  1020,  which  provides  an  upper  bound  for  a  , 

1  1 

so  that  the  optimal  shapes  are  only  slightly  curved  over  S'OS. 

This  one  relation,  a  =  f (if),  is  all  that  is  needed  to  complete  the 
solution  since  the  factor  A  drops  out  of  the  expressions  for  the  drag  co¬ 
efficient  and  the  ratio  of  arc  length  to  chord.  Thus,  by  (36), 

CD  =  CD(a,K)  =  D*/yo  ,  (57) 

and,  denoting  the  ratio  of  arc  length  to  chord  by  k, 

k  -  k(a,K) r s0/y0  .  (58) 

The  evaluation  of  (57)  and  (58)  (in  which  D  ,  yQ,  and  so  are  given  by  (52), 
(53),  and  (54)  )  for  a^  =  £{k)  gives  a  parametric  representation  of 
versus  k.  This  is  plotted  in  Fig.  7,  where  Cn  =  2n7(ir+4)  is  the  drag 

coefficient  of  a  flat  plate  in  infinity  cavity  flow  (e.g. ,  see  Lamb(1932)  ). 

At  k  -*■ oo.  it  can  be  shown  from  (57)  and  (58)  that  C^^4+’n)/(ir+8k). 

Minimum  drag  profiles  for  various  values  of  k  are  shown  in  Fig.  8 . 
These  profiles  are  found  by  numerically  integrating  dz  in  (17)  and  are 
quite  similar  to  Lavrentieff's  profiles  discussed  earlier  in  Section  1;  how¬ 
ever,  by  expanding  yQ  in  (53)  and  s^  in  (54)  for  small  a ,  it  can  be 
shown  that  for  (k-1 )  «  I , 
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cd  =  tet {1 ' Y{k_1 }*  +  °(k"1 )}  • 

where 

X 

Y  =  4{2(ir+4)}"2  ~  1.0584 

for  the  Lavrentieff  profiles  (a^  =  0),  and 

Y  =  4{(9ir+64)/4(ir+4)(3T!+16)P  ~  1.1641 

for  the  profiles  in  Fig.  7.  Thus,  for  k  close  to  unity,  the  drag  coefficients 
of  the  profiles  in  Fig.  7  are  slightly  less  than  those  for  Lavrentieff's  profiles. 

The  cases  N  =  2,  3,  .  .  .,  could,  in  principles,  be  carried  out 
along  similar  lines  and  should  result  in  improved  drag  coefficients  for  a 
given  k  =  s^/y^.  The  numerical  examples  of  Part  I,  in  which  the  exact 
solutions  to  the  variational  problem  are  known,  indicate  that  expansion  in 
Fourier  series  is  a  very  effective  method,  at  least  for  the  case  of  quadratic 
functionals.  Whether  the  same  holds  true  for  the  present  problem,  in  which 
the  functional  is  of  a  different  type,  remains  to  be  seen. 
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